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In this thesis we explore the character of finite-time singularities that are possible to arise 
in FRW universes and in specific braneworld configurations. 

We begin with a review of the structure and properties of RW-spaces. We then give 
an adaptation of a completeness theorem and analyse complete model universes from the 
recent literature. 

We continue with the construction of a classification plan of singularities that may 
appear in FRW universes. Our tool for this quest is the contrapositive of the complete- 
ness theorem which translates in necessary conditions for singularities. The classification 
derived in this way is based on the behaviour of the Hubble rate and offers a first insight 
into the nature of singularities. We further extend and complete this classification with 
the use of the Bel-Robinson energy. We also associate the existence of a closed trapped 
surface with a divergence in the Bel-Robinson energy and use this relation to determine 
the evolution of various specific cosmologies. 

Next we study models from the recent literature that offer the appropriate physical 
conditions for the appearance of singularities predicted by our classification. This study 
shows that our classification not only accommodates all known types of singularities but 
it also predicts new ones. 

Finally, we focus on braneworld models consisting of a three-brane embedded in a five- 
dimensional bulk space that is inhabited by a scalar field or a perfect fluid. Our goal is to 
study the singularities that may occur within a finite distance from the brane. Our analysis 
that is performed with the method of asymptotic splittings shows that these braneworlds 
can exhibit three main types of finite-distance singularities. 
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Chapter 1 
Introduction 



Ever since the fundamental discovery of the big bang and big crunch singularities in the 
general relativistic FRW universes of dust and radiation a lot of effort has been focused on 
finding ways to avoid these philosophically and physically puzzling features of spacetime 
structure. Indeed, spacetime singularities constitute a challenging issue philosophically 
since they imply that the immense universe emerged a finite time ago through a region of 
spacetime of infinite curvature and zero size, as well as physically since general relativity 
fails to describe such a state of spacetime. A common way out of this situation is to modify 
general relativity to obtain theories of gravity that do not possibly contain singularities. A 
different attitude is to accept that singularities are essential features of our universe and 
investigate further conditions that trigger their existence as well as those that describe 
their nature. 

The singularity theorems of general relativity [1] gave appropriate conditions that led to 
the existence of singularities in generic spacetimes; they did not, however, provide detailed 
information about their nature. Uncovering the nature of spacetime singularities remains 
therefore an open issue in cosmology and this thesis focuses on the precise pathologies that 
may determine the character of singularities in the class of isotropic spacetimes. Before we 
analyse more carefully our goals, let us briefly review the development of recent research 
on the issue of finite-time singularities that partly motivates our own research reported 
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herein. 

An investigation of the recent literature actually shows that interest in the issue of 
spacetime singularities has been renewed through the construction of various cosmological 
models with finite-time singularities. Among these are those similar to the standard big 
bang and big crunch singularities ones albeit appearing in universes containing different and 
even 'exotic' forms of matter, as well as other kinds not met in the standard cosmological 
models studied in the usual textbook expositions of cosmology. 

Interestingly, these novel types of singularities sometimes involve the presence of quite 
an exotic type of matter (cf. [2]- [17]), a characteristic example of such a singularity being 
the so-called big-rip singularity (cf. [3], [H], [SB], [IB], 03], [SI], [SI], [52], P])- The 
study of some of these newer types of singularities started off as soon as data from the 
WMAP (cf. [51], [55]) became available. These data indicate that the universe is currently 
undergoing an accelerating phase typically attributed to a form of dark energy that may be 
described by a linear equation of state with a parameter w converging, however, towards 
the value w = —1. The case with w < — 1 is not ruled out by the data [3J and this provides 
a surprising implication for the future of the universe which had not been appreciated 
earlier: once such a dark energy component becomes dominant, the universe is destined to 
face up with a violent end towards which the scale factor, the Hubble expansion rate, the 
horizon distance, the density and the pressure all experiencing a divergence. In this sense 
every bound structure will be ripped apart. 

The type of dark energy that leads to the emergence of the big-rip type of singularity 
mentioned above has some awkward characteristics (see Section 4.3) but it nonetheless 
sparkled new interest on finite-time singularities. In particular, research focused (cf. [4"9] . 
|50j, [56j-[66j) in finding paradigms of cosmologies that exhibited future finite-time sin- 
gularities potentially distinct from big crunches and rips but still feasible to arise in an 
expanding phase of a model universe which did not require exotic forms of matter. This 
led to the discovery of sudden and other milder types of singularities. These singularities 
are characterised by the fact that the second or higher order derivatives of the scale factor 
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diverge while the scale factor itself and the Hubble rate remain finite. 

The models described in the preceding paragraphs played the important role of pro- 
viding the appropriate physical circumstances for the appearance of such singularities. 
But still they provide mere examples not general theorems that could capture the unify- 
ing reasons for the occurrence of such spacetime singularities. In fact, the study of an 
increasing number of models accommodating finite-time singularities leaves one with the 
impression that different singular behaviours are totally determined by (or connected to) 
the particularities of each model. 

This situation, however, provides us with a fertile ground for the formulation of some 
basic questions such as: What do all these different singular cosmologies have in common? 
Is there some underline reason that unifies all these 'phenomenologically' different singular 
behaviours? Are there other types of finite-time singularities which are possible in these 
models but have not yet been uncovered? What are the general criteria that predict and 
describe the possible kinds and nature of finite-time singularities? 

This thesis aims to provide answers to such questions. We will show that singular 
isotropic universes are characterised by the basic property that their Hubble expansion rate 
fails in one way or another to be an integrable function of the proper time for an infinite 
time-interval. The various ways in which this non-integrability can be implemented leads 
to a classification of finite-time singularities that consists of three possible types which 
accommodate all the recently introduced singularities mentioned above [HZ], [HE]- This 
result offers first evidence for a more detailed description of the nature of singularities in 
isotropic spacetimes. 

Note that such a 'classification' of singularities is totally effected by the geometry of 
the spacetime in question; in no way do the various matter fields influence its construction. 
A more subtle examination of the character of a singularity implies however that matter 
fields may strongly affect the asymptotic behaviour of these model universes and hence 
play a role in changing the nature of their singularities. This type of classification needs 
therefore to be refined by somehow incorporating the various matter fields into the picture. 
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This refinement results in a more detailed classification scheme that sharpens our view of 
the nature of singularities [69j-[71j. 

Our approach in constructing a classification of singularities is a geometric one. In 
particular, for the Hubble rate classification mentioned above we make use of the contra- 
positive of the completeness theorem of [72], [73], while for the second we use the notion of 
the Bel-Robinson energy (an invariant geometric quantity representing the slice-projected 
energy of the gravitational field). This latter approach makes our classification conducive 
to more general spacetimes as well as more flexible since now more functions describe the 
possible different types of singularities and trace the properties and structure of spacetime. 

After the completion of the classification of singularities in isotropic universes, this the- 
sis focuses on the analysis of singularities that occur in the different context of braneworld 
universes. The study of these models (originally motivated by string theory) offers many 
interesting scenarios in cosmology and particle physics (see Section 5.1 for a short review). 
Our interest in this recently issued research area originates from the fact that these models 
exhibit finite- distance singularities that have many analogies with the finite-time singular- 
ities predicted in our classification and not met in simpler, more standard models. 

For instance, the specific braneworld model of [71] captured our interest precisely 
because it involves the appearance of a finite-distance singularity. This model consists of a 
flat three-brane sitting at a fixed position of an extra spatial dimension, and a scalar field 
propagating throughout the whole of the enveloping five-dimensional space. At the vicinity 
of the singularity occurring within a finite-distance from the brane, the energy density of 
the scalar field becomes infinite and this may offer an explanation to the small, observed 
value of the vacuum energy (actualised in this model by the scalar field). 

We further explore whether this singularity continues to appear when we consider a 
curved instead of a flat brane. Our preliminary results indicate that the singularity is in 
fact removed at an infinite distance from the brane [75]. This work is further extended 
here by examining other types of finite-distance singularities that become possible once we 
consider a different form of bulk matter such as a perfect fluid. 
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By the use of the method of asymptotic splittings developed in [76] , we find that there 
are three possible types of finite-distance singularities, triggering an infinite leak of energy 
of the matter that inhabits the background space. 

A more detailed plan of this thesis is as follows: 

In Chapter 2, after a review of the basic properties of Robertson- Walker spaces we 
give an adaptation of the proof of the general completeness theorem of [72] for the specific 
case of such spaces and analyse complete model universes from the recent literature. 

In Chapter 3, we use the contrapositive of the completeness theorem [72] , [73] to derive 
necessary conditions for singularities in isotropic spacetimes. These in turn provide us with 
a first classification plan of singularities. This plan is later extended and completed with 
the use of the Bel-Robinson energy. We also associate the existence of a closed trapped 
surface with a divergence in the Bel-Robinson energy and use this relation to determine 
the evolution of various specific isotropic models. 

In Chapter 4, we exploit model universes appearing in the recent literature to demon- 
strate how some particular members of our classification indeed arise during the evolution 
of such models. It follows from the analysis in this chapter that our classification scheme 
accommodates all known types of singularities. In fact, an examination of the literature 
shows that many new types of singularities -predicted by our classification- have not yet 
been embodied in appropriate cosmological models. 

In Chapter 5, we present an asymptotic analysis of braneworld models that consist 
of a three-brane embedded in a five- dimensional bulk space filled with a scalar field or a 
perfect fluid. Our analysis, performed with the method of asymptotic splittings, shows 
that these models are constrained to exhibit only three main types of singularities within 
a finite distance from the position of the brane which is taken to be at the origin. 

In Chapter 6 we conclude and discuss projects of future work. This thesis has also 
five appendices. In Appendix A, we explain the notation and conventions that we use 
throughout, and in the appendices that follow we give some background of some of the basic 
concepts that we use in this work such as the extrinsic curvature and the Gauss-Codazzi 
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formalism (Appendix B), the Bel- Robinson energy (Appendix C), the energy- momentum 
tensor of a scalar field (Appendix D) and the method of asymptotic splittings (Appendix 
E). 

We declare that all statements in this work correspond to original research conducted 
by the author except in all those places where it is otherwise stated and proper citations 
are then given. 



Chapter 2 



Completeness and Robertson- Walker 
geometry 

In this Chapter we explore the character of Robertson- Walker spaces, the basic spacetimes 
used and studied in this work. We start by reviewing, in Section 2.1, the motivation 
behind the adaptation of such spaces, their basic properties and structure as well as some 
simple features of singularities that they may exhibit. In Section 2.2, we first review the 
assumptions of the completeness theorem stated and proved in [72] for generic spacetimes 
and then give a version of that proof for the case of a Robertson- Walker space. The key 
idea of the proof of that theorem is to study the geodesic equations satisfied by every causal 
geodesic; this is a different approach from the one taken in the proof of the 'antipodal' 
singularity theorems, where the idea is to study instead the Raychaudhouri equation which 
describes the expansion of a bundle of geodesies. The completeness theorem applied to 
the case of a Robertson- Walker space guarantees that the integrability of the norm of the 
extrinsic curvature, or equivalently of the Hubble expansion rate, on an infinite-proper 
time interval suffices for the completeness of such spaces. Finally, in Section 2.3 we analyse 
complete universes from the recent literature. 
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8 CHAPTER 2. COMPLETENESS AND ROBERTSON-WALKER GEOMETRY 

2.1 Structure and properties of Robertson- Walker spaces 



Since observations on the large-scale structure of the universe are limited, the construction 
of cosmological models includes at a most basic level philosophical ideas that we accept 
as principles. The cosmological principle plays a central role in the construction of model 
universes and can be viewed as a generalisation of the Copernican principle according to 
which the earth does not occupy a privileged position in the solar system. The cosmological 
principle assumes that in turn the solar system, our galaxy, and our local group of galaxies 
do not occupy a privileged position in the universe. 

The cosmological principle therefore implies that the universe on a large enough scale 
ought to be spatially homogeneous. A spacetime is spatially homogeneous if at any instant 
of time, represented by a spacelike hypersurface St, any two points on S t are equivalent. 
More precisely, a spatially homogeneous spacetime can be foliated by a one-parameter 
family of spacelike hypersurfaces S 4 such that for each t and for any pair of points p, q on 
S t , there exists an isometry that maps p into q (cf. [77], pp. 92-93). This implies that the 
curvature of S 4 cannot vary from point to point i.e., the spacelike slices S t are surfaces of 
constant curvature. 

Apart from imposing homogeneity, the cosmological principle also requires isotropy 
about every point in spacetime; this means that there are no privileged directions, or 
equivalently, that the universe is spherically symmetric about every point. 

The strongest evidence for isotropy is provided by the discovery of Penzias and Wilson 
[78] that the universe is filled with a microwave background radiation of about 3 K that 
is highly isotropic. Since isotropy implies homogeneity [79], this fundamental discovery 
provided evidence for the validity of the cosmological principle. 

The requirement of isotropy further implies that there exists a class of privileged ob- 
servers, the isotropic observers, that move along a congruence of timelike geodesies with 
tangent vectors orthogonal to the homogeneous slices S t . These timelike geodesies do not 
intersect, except possibly at singular points, so that through every non-singular point in 
spacetime passes at most one isotropic observer. The orthogonality between these geodesies 
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and the slices E f is a necessary feature of isotropic spacetimes: if these geodesies did not 
cut the surfaces H t orthogonally, then we would be able to construct a preferred spatial 
direction in violation of isotropy (cf. [77J, p. 93). 

The suitable metric for the description of an isotropic spacetime is the Robertson- 
Walker metric given by 

ds 2 = -dt 2 + a 2 (t)dtt 2 , (2.1.1) 

where dQ 2 is the metric of the three-spaces of constant negative, zero or positive curvature. 
The non-zero constant curvature can be normalised to be ±1, (so that dQ 2 is the metric of 
the unit three-sphere or of the unit hyperbolic three-space respectively, and the trichotomy 
of possibilities reads as follows: 

ds 2 = -dt 2 + a 2 (t)(dr 2 + f 2 (r)(d6 2 + sin 2 6d(j) 2 )), (2.1.2) 

where a(t) is the scale factor, and 



r, k = 



smr 



k=l (2.1.3) 



sinhr, k = —1, 



with r G [0, oo) if k = 0,-1, whereas r G [0, 2ir] if k = +1. We call a universe with 
constant curvature k equal to —1,0, or +1, an open, flat, or closed universe respectively. 
Note that the spacelike slices of a flat or open universe are of an infinite extent in contrast 
to those of a closed universe which are finite. 

As we have already remarked already in a isotropic universe we are always able to 
construct a privilege class of isotropic observers. Each of these observers can be thought 
of as a particle of a cosmic fluid filling spacetime; the world lines of the observers become 
then the flow lines of the fluid. Because of the orthogonality between these flow lines 
and the spacelike slices E t , the spacelike coordinates of a particle in the fluid are constant 
-co-moving coordinates- along these lines. This cosmic fluid models the matter content of 
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the universe and the scale factor describes the separation of neighboring flow lines (cf. pQ 
p. 136). In the simplest case the energy-momentum tensor of the cosmic fluid is thus taken 
to be that of a perfect fluid with density \i and pressure p which are functions only of the 
(proper) time. The resulting Einstein-Friedman equations can be written in the form 

3H 2 = /i + 3-^, (2.1.4) 

cr 

- = --(// + 3p), (2.1.5) 
a 6 

while the conservation of the energy-momentum tensor gives 



{i + 3H(n + p) = 0, (2.1.6) 

where H = a/a is the Hubble rate. 

In the absence of exotic forms of matter, we assume that \x > and p > for all 
t. Then Eq. (12.1.51) shows that a cannot be constant and thus the universe is either 
expanding or contracting. The discovery of the expansion of the universe by Hubble in 
1929 offered a very interesting implication for its history: because its density decreases as 
the universe expands (see Eq. (I2.1.6P ). it follows that it should have been infinite in the 
past as a(t) — ► 0; the universe must have therefore began its existence with a singularity. 
In fact, the singularity theorems predict that this is inevitable in all models satisfying 
fi + 3p > 0. 

The first cosmological models with a RW geometry that were found to possess such a 
singularity were the dust and radiation-filled universes. Dust is a pressureless fluid that 
provides a good approximation for a matter component comprised solely by non-interacting 
galaxies. This follows because the observed relative velocities of galaxies in groups are small 
and hence galaxies are interpreted as particles of a pressureless fluid. On the other hand, an 
adequate approximation to the matter content of the universe in an early epoch is provided 
by radiation which is described by taking p = /i/3. Away from an initial singularity open 
and flat universes typically expand forever, whereas, closed ones end their existence with 
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another singularity that arises within a finite time from the initial big bang. 

2.2 Completeness theorems 

Although we are going to focus on exclusively isotropic spacetimes, we refer briefly in 
this Section to the basic properties of general spacetimes used in the formulation of the 
completeness theorem of [72]. In the end of this Section, we give an adapted version of 
the original theorem and proof of [72] , for our case of interest -isotropic spacetimes. We 
will refer to the original version of the completeness theorem ( Theorem 12.2.11 below) later 
in Chapter 3 when we shall discuss the possible extensions of our work to more general 
spacetimes. 

We begin by giving some basic definitions. Consider a spacetime (V, g) with V = 
M. xX, X = (to, oo), where M. is a smooth manifold of dimension n and ( n+1 )g a Lorentzian 
metric which in the usual n + 1 splitting, reads 

( n+ Vg = ~N 2 dt 2 + g l] (dx i + (3 l dt){dx j + (3 j dt) (2.2.1) 

Here N = N(t, x l ) is called the lapse function, (3 l (t, x^) is called the shift function and the 
spatial slices M. t (= M. x {t}) are spacelike submanifolds endowed with the time-dependent 
spatial metric g t = g^dx % dx 3 ' . We call such a spacetime a sliced space [SU]. A sliced space 
is time-oriented by increasing t and we choose X = (to, oo) (or X = (— oo, to)) when we 
study the future (past) singularity behaviour of an expanding (contracting) universe with 
a singularity in the past (future), for instance at t = < to (t = > t ). However, since t 
is just a coordinate, our study could apply as well to any interval IcR. 

A natural causal assumption for (V, g) is that it is globally hyperbolic. This implies the 
existence of a time function on (V, g). In a globally hyperbolic space, spacetime splits as 
above with each spacelike slice Ait being a Cauchy surface; this means that each timelike 
and null curve without end points cuts M. t exactly once [81]. In such a spacetime therefore 
the future and history of the universe can be predicted or retrodicted from conditions at 
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the instant of time represented by Ait (cf. [77], pp. 200-202). 

Following the definitions given in [72] , we say that a sliced space has uniformly bounded 
lapse if the lapse function N is bounded below and above by positive numbers N m and 

N M , 

< N m < N < N M - (2.2.2) 

A sliced space has uniformly bounded shift if the g% (t G X) norm of the shift vector (3, 
projection on the tangent space to Ait of the tangent to the lines {x} x X, is uniformly 
bounded by a number B. 

A sliced space has uniformly bounded spatial metric if the time-dependent metric g t = 
gijdx l dxi is uniformly bounded below for all t G X by a metric 7 = g to , that is there exists 
a number A > such that for all tangent vectors v to Ai it holds that 

A lij v i v j < gijvV. (2.2.3) 

Definition 2.2.1 A sliced space (V,g) with uniformly bounded lapse, shift and spatial 
metric is called regularly sliced. 

Denoting by ViV the space gradient of the lapse N, by the extrinsic curvature of 
Ait, and by \K\ 2 g the product g am g^K^K^^ we have the following theorem of [72] which 
gives sufficient conditions for geodesic completeness: 

Theorem 2.2.1 Let (V, g) be a sliced space such that the following assumptions hold: 
CI (V,g) is globally hyperbolic 
C2 (V, g) is regularly sliced 

C3 For each finite t\, the space gradient of the lapse, \ \7N\ g , is bounded by a function of 
t which is integrable on [£i,+oo) 

C4 For each finite t±, \K\ g is bounded by a function oft which is integrable on [ti,+oo). 
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Then (V, g) is future timelike and null geodesically complete. 

It was later proved in [80] that in a regularly sliced spacetime, condition CI is in fact 
equivalent to the condition that each slice of (V, g) is a complete Riemannian manifold. 

Turning back now to isotropic spacetimes we see from the definition of a sliced space 
and the form of the RW metric (Eqs. (12.1.21) . (12.1.31) ). that an isotropic spacetime is a 
sliced space with N — 1, (3 — 0. Following the Definition (12.2.11) we note that: 



Remark 2.2.1 An isotropic spacetime is regularly sliced iff the scale factor is a bounded 
from below function of time. 

We can now reformulate the completeness theorem 12.2.11 so as to apply directly to the 
case of an isotropic spacetime. 

Theorem 2.2.2 (Completeness of RW spacetimes) Every globally hyperbolic, regu- 
larly sliced RW space such that for each finite t\ the norm of the extrinsic curvature is 
integrable on [t\, +oo) ; is future timelike and null geodesically complete. 

Proof. 

Consider a causal geodesic in an isotropic universe. The tangent vector to these 
geodesies u a = dx a /ds, where x a = (x°, x 1 , x 2 , x 3 ) = (t,r,9,4>) and s is the proper time or 
an affine parameter measured along a timelike or null geodesic respectively, satisfies the 
geodesic equations 

_____ _i_ p a ^_^±_ _ n (2 2 41 

ds z ds ds 

The a = component of these equations for this type of universe is 

d 2 t _ n dx % dx j 

^ + r °TO = °- < 2 - 2 ' 5 > 

where i, j = 1, 2, 3 (the rest of the Christoffel symbols and Tq vanish). 
Dividing Eq. fT2T2~5|) by (dt/dsf we find, 

ds d ( dt\ Q dx^xi = 

dtdt \ds) ij dt dt { ' ' ' 
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Setting y = dt/ds and denoting by y the derivative of y with respect to t, the above 
equation is written in the form 

(2.2.7) 



y 

Integrating over the interval [t\,t], we find 

y(t) 



y o ax ax 

ij ~dT~dT 



In 



y{ti 



(2.2.8) 



where we have set v l = dx % jdt. Since r°,- = i^y □ we can write Eq. (I2.2.8|) as follows: 







y(ti) J h 3 



(2.2.9) 



Substitution then in Eq. (3.16) of [72] gives 

In yjfi < f\ K \ dt = ^ f 

y{ti) J tl J tl 



f 




a 


hi 


a 



dt, 



(2.2.10) 



where in the last equation we have used Eq. (IB. 1. 15[) from Appendix B. Using the hypoth- 
esis of the integrability of \K\ g , we see that y = dt/ds is uniformly bounded. This means 
that ds/dt is bounded away from zero and hence the length of the geodesic given by 



ti 



ds 
~dt 



dt 



(2.2.11) 



is infinite. This implies that the spacetime is causally geodesically complete. ■ 

Past completeness is derived by exchanging future with past and \pi, +oo) with (—00, ti]. 
It follows from the above theorem that a globally hyperbolic and regularly sliced isotropic 
spacetime is complete provided that the spacelike slices H t do not curve too much as they 
are seen by an observer located outside T, t in spacetime, for all t G [t±, +00). 

Now since the norm of the extrinsic curvature in an isotropic spacetime is proportional 



1 Using Eqs. (|B.1.14|) - (|B.1.14|) we can verify that r? x = da = K n , T° 22 = hap = K 22 , ^33 = 
aaf 2 sin 2 9 — K33. 

2 As it was shown in [72] a future-directed causal curve C : s 1— * C(s) can be reparametrised by t: if 
n = (— 1, 0) is the timelike normal to Mt then g (^f, n) — — lj < and hence 4| > 0. 
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to the Hubble expansion rate H(t) (see Appendix B, Eq. ( IB. 1.151) ). we conclude that such 
a spacetime is complete if it exhibits an expansion rate that is an integrable function of 
proper time on \b\, +oo). 

In the following Section, we will see that there exist complete FRW models that have 
an expansion rate which is non-integrable on \t\, +oo). This implies that the integrability 
of the Hubble rate on an infinite proper-time interval is a sufficient but not a necessary 
condition for completeness. In order to prove completeness of these universes we take a 
different approach in which we make use of the Bel-Robinson energy. 

As we discuss in Appendix C, the Bel-Robinson energy is a kind of energy of the 
gravitational field projected in a sense to a slice in spacetime. In particular, for a RW 
metric the Bel- Robinson energy at time t, £>(£), is given by 



where C is the constant volume of (or in, in the case of a non-compact space) the three 
dimensional slice at time t and \E\, \D\ are the norms of the two electric tensors, namely, 



It can be proved that a closed or flat, expanding at FRW universe that has a 
bounded Bel- Robinson energy is causally geodesically complete jH2]. Therefore complete 
FRW models which fail to have a Hubble expansion rate that is integrable on an infinite 
proper-time interval may still manage to be complete provided that they have a bounded 
Bel-Robinson energy. In the following Section we will illustrate this result with the analysis 
of various models. 




(2.2.12) 




(2.2.13) 



(2.2.14) 
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2.3 Complete FRW cosmologies 

Completeness is not a property only shared by 'normal' cosmologies. Although later in 
Section 4.3, we meet isotropic universes filled with dark energy that end their existence in 
a big-rip singularity, this need not always be the case. There are many examples of such 
'exotic' universes containing dark energy that still manage to evade the big-rip singularity 
and exist forever. This is achieved either by assuming that the dark energy component 
satisfies simultaneously two different equations of state, or that it satisfies a generalised 
equation of state that includes matter ranging from Chaplygin gasses to perfect fluids. We 
shall presently analyse examples of complete universes belonging to both of these categories 
in turn. 

For the first case, consider the flat FRW universe studied in [83] (see also [84]). In this 
case the dark energy component satisfies a 'phantom' equation of state of the form 

p = wfi, with w = w(t) and w < — 1, 

as well as an equation of state of a Chaplygin gas (cf. [85], [86] ) 

p = , where A > 0, constant, 

(similar results will hold when w G (—1,-1/3) -/c-essence models, see, e.g., [El]). Sub- 
stituting only the latter equation of state in the continuity equation, Eq. ( 12. 1.6ft . and 
integrating we find 

^{t)=A + (j4-A)(^j . (2.3.1) 

We can now take into account the phantom equation of state and reform the above equation. 
This is done by calculating, with the use of both equations of state, the time-dependent 
parameter w(t), 

w(t) = -4- (2.3.2) 
Substituting this back in Eq. (12.3. If) . we find the following form for the energy density, 
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fi(t) = Vo 2 [-w + (1 + w )(a /a(t)) 6 f 2 , (2.3.3) 

where A = — u; /i 2 and Wq < — 1, which we can now insert in the Friedmann equation, Eq. 
(12.1.41) . and derive the following solution for the scale factor: 

a(t) = (Ce Cl( '-' o) ) 1/6 , (2.3.4) 

where C and C\ are positive constants. This model admits cosmic acceleration, a > 0, 
and a(t) — > oo as t — > oo and it is geodesically complete. The Hubble parameter takes the 
constant value, 

H = ^f. (2.3.5) 
o 

Here also the Bel-Robinson energy remains always finite and hence by the completeness 
theorem in |82j we conclude that the model is geodesically complete. 

As we mentioned in the beginning of this section, the second technique to trigger the 
avoidance of a big-rip singularity in models containing dark energy is to assume that the 
dark energy satisfies a generalised Chaplygin equation of state. This was done in [81] for 
a flat universe. In particular, the equation of state is here chosen so that 

p = -fT a [C + (fi 1+a - C) a/{1+a) ], (2.3.6) 

where C = A/(l + w) — 1, 1 + a = 1/(1 + w), and a is a real parameter. Note that in the 
limit w — > the gas reduces to a Chaplygin gas satisfying p = —An~ a , whereas as A — > 
it takes the form of a perfect fluid with equation of state p = wfi. The scale factor is given 
by the form 

a(t) = (de- ^-^ + C 2 e c ^- to) ) 2l \ (2.3.7) 

where Ci, C*2 and C3 are parameters depending on ao (the initial value of the scale factor), 
A, w and B, with B being a positive integration constant. At the asymptotic limits t — > to 
and t — > 00, we find that H tends to suitable constants, that is it remains finite on [to, 00). 
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The complete character of this model is guaranteed by a clearly bounded Bel-Robinson 
energy. 

Assuming the existence of more conventional matter contents such as radiation and a 
cosmological constant in a closed FRW universe, we obtain the following solution originally 
given in [87], 

r r- I 1 / 2 

a{t) = a % 1 + exp(V2t/ai) , (2.3.8) 

where aj is a constant that bounds a(t) from below making this universe regularly hyper- 
bolic. One then finds that H becomes zero as t — > — oo, it is integrable on (— oo, 0] and 
the Bel-Robinson energy is bounded. Thus we conclude that this universe is geodesically 
complete, either using the result of [82] or alternatively by the completeness theorem 12.2.21 



Chapter 3 



Classification of singularities 

3.1 Introduction 

Up to now we studied situations where geodesic completeness plays a major role. From now 
on we shall focus on the other important issue, namely, that of gravitational singularities. 
Spacetime singularities characterise solutions describing cosmological and gravitational col- 
lapse situations in general relativity and in other metric theories of gravity. Understanding 
the various aspects of spacetime singularities constitutes the so-called singularity prob- 
lem. Viewed from a mathematical perspective the singularity problem apart from being 
appealing because of its fundamental nature, it is also a pressing one, ever-present in man- 
ifolds and metrics of arbitrary form under general topological and geometric assumptions 
(see Section 3.4). From a physical perspective the singularity problem attracts attention 
because of the wide variety of cosmological models in general relativity and other metric 
theories of gravity that exhibit singularities that need to be analysed and understood (see 
Sections 4.2-4.5 for a detailed analysis of singular cosmological models in general relativity). 

Even though spacetime singularities represent an extensively studied area from both a 
mathematical as well as a physical point of view, there are still many important unresolved 
issues in this general area. It is essential to find appropriate tools that will serve to define 
and analyse the structure and nature of spacetime singularities. This will hopefully lead to 
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a more basic understanding of all those qualities that make up the character of singularities 
along with a more detailed classification of them. But there exists a risk in performing such 
a project: mapping out only a part of the field comprised of all possible types of singularities 
would provide us with misleading conclusions about the generality of the types found since 
we would then be shedding light only on some aspects of the field while others would 
remain concealed [88]. It is thus important to attempt to plot a more or less complete 
classification scheme. Any such scheme ought to acknowledge both the mathematical and 
the physical qualities of a spacetime singularity. Among the mathematical qualities the 
most important one is that of geodesic completeness/incompleteness and it should therefore 
play a key role in the planning of a classification. On the other hand, physical qualities 
describe the way matter fields behave on approach to the singularity and are expected to 
contribute significantly to our understanding of its nature. 

We begin the analysis of singularities in this chapter by focusing on those types of 
singularity that may arise in the simplest class of isotropic spacetimes. The tools that we 
use are of a geometric nature. Our geometric approach on one hand makes the resulting 
classification scheme easily adjustable to a further future study of more complicated space- 
times, and on the other it helps us find complementary results and generic properties of 
the class of spacetimes being analysed. 

Our attack to a realisation of a complete classification scheme as discussed above 
is implemented in two stages. In the first stage, we exploit the contrapositive of the 
completeness theorem Theorem 12.2.11 in the effort to derive necessary conditions for the 
existence of singularities in the case of isotropic cosmological models. This provides us 
with a classification of singularities that is solely based on the behaviour of the Hubble 
expansion rate H. 

For various reasons to be given later, the classification depicted in the first stage does 
not provide us with all details necessary for a complete classification scheme referred to 
above. We therefore need a complementary tool in order to bring to surface all the other 
missing factors and complete our plan. 
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This leads us to the second stage which consists of a classification scheme based on the 
introduction of a new invariant, the Bel- Robinson energy, a geometric quantity built out 
of combinations of various curvature components that traces successfully all those required 
qualities missing from the first classification. 

Apart from being involved in the second stage of our classification scheme, we show 
that the use of the Bel-Robinson energy in combination with the notion of a closed trapped 
surface induces important information about the dynamical evolution of isotropic models. 
In particular, we prove that the existence of a closed trapped surface can be traced in the 
dynamical behaviour of the Bel-Robinson energy. We show by analysing various examples 
from the literature that when the existence of a closed trapped surface leads to a singularity 
the Bel-Robinson energy also diverges, whereas in the opposite case the Bel-Robinson 
energy remains finite and signatures the complete character of the model. 

The plan of this chapter is as follows: In Section 3.2, we classify singularities according 
to the behaviour of the Hubble rate. In Section 3.3 we analyse the necessity for refining 
our classification scheme: In Subsection 3.3.1 we derive a complete classification scheme 
based on the possible behaviours of three functions the Hubble rate, the scale factor and 
the Bel-Robinson energy. In Section 3.4 we analyse the interplay of closed trapped surfaces 
and the dynamical behaviour of the Bel- Robinson energy: In the Subsection 3.4.1 we show 
how the existence of a closed trapped surface is linked to a possible divergence in the 
Bel-Robinson energy and in Subsection 3.4.2 we illustrate this result for various model 
universes. Finally, in Section 3.5 we discuss our conclusions. 

3.2 Hubble parameter singularities 

Theorem 12.2.21 provides us with sufficient conditions for completeness and we can therefore 
use its contrapositive to derive necessary conditions for the existence of singularities. This 
will also give us an insight into the character of possible singularities. In particular, consider 
a singular, globally hyperbolic and regularly sliced isotropic universe. Then according to 
Theorem 12. 2. 21 there is a finite time t\ for which H fails to be integrable on the proper time 
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interval [t\, oo). In turn, this non-integrability of the expansion rate H can be implemented 
in different ways and we arrive at the following result which provides first evidence for the 
nature of finite-time singularities, based entirely on the behaviour of the Hubble parameter 

Theorem 3.2.1 Necessary conditions for the existence of finite time singularities in glob- 
ally hyperbolic, regularly sliced isotropic universes are: 

51 For each finite t, H is non-integrable on [t\,t], or 

52 H blows up in a finite future time, or 

53 H is defined and integrable for only a finite proper time interval. 

Condition Si describes a big bang type of singularity when H blows up at t\ since then it 
is not integrable on any interval of the form [ti,t], t > t\ (the condition of regular slicing 
is violated in this case but the scale factor is bounded from above). However, under S\ we 
can have other types of singularities: Since H{r) is integrable on an interval [ti,t], if H{t) 
is defined on [t\,t], continuous on (ti,t) and the limits lim T _^+ H{r) and lim r _,. t - H(t) 
exist, the violation of any of these conditions leads to a singularity that is not of the big 
bang type discussed previously. 

Condition S2 describes a future singularity at t s (t s > t±) characterised by the diver- 
gence of H at t s . Condition S3 may lead to a future singularity at t s with H being finite 
on [t a ,t s ] where t a > t 1 but for this to be a genuine type of singularity, in the sense of 
geodesic incompleteness, one needs to demonstrate that the metric is non-extendible to a 
larger interval. 

Note that these three conditions are not overlapping: clearly S± and S2 are not over- 
lapping with S 3 but also Si is not overlapping with S 2 - Indeed Si is not implied by S 2 
for if H blows up at some finite time t s after ti, then it may still be integrable on [ti,t], 
ti < t < t s . 

There are many examples of singularities met in the literature that belong to the 
types predicted by the above theorem; they are analyzed in Chapter 4 in detail. They 
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usually describe flat isotropic universes with various components of matter. For example, 
the phantom dark energy models of [51] exhibit a future singularity which falls in the S 2 
category. The standard big bang singularities as well as the sudden singularity of (56] fall 
in the Si category. Other sudden singularities (which do not have a blow up of H at t = 
(cf. [52J)) and the inflation model of [89] both fall in the S3 category (for a detailed analysis 
of these models see Chapter 4). 



3.3 The (S,N) pairs 

Although the classification based on the Hubble parameter is a first broad step to clearly 
distinguish between the various types of singularities that can occur in isotropic universes, 
it does not bring out some of the essential features of the dynamics that may differ from 
singularity to singularity. To see this we note that S 3 encloses altogether dissimilar types of 
singularities each characterised by a distinct pathology, for example in higher order deriva- 
tives of H, and as we have already mentioned we also accommodate distinct singularity 
types under Si. In addition, condition S 2 includes both a collapse singularity, where the 
scale factor a —>■ as t —>■ t s , and a singularity with a —>■ 00 as t —>■ t s . This type of 
degeneration is unwanted since it impedes the construction of a complete classification of 
singularities. It is therefore necessary to extend and refine this classification. 

We begin by including in our classification the behavior of the scale factor. That it is 
only necessary to include this behaviour is seen most clearly by noticing that only in this 
way one may consistently distinguish between initial and final singularities, or between past 
and future ones, as it has been repeatedly emphasised by Penrose in his Weyl curvature 
hypothesis (see for example [9U] PP- 765-769). According to this physical conjecture, the 
Weyl curvature which acts as a distorting source for matter is constrained to be almost 
zero at initial singularities. In contrast, at final singularities it diverges to infinity. 

Taking into account the behaviour of the scale factor in our singularity classification 
scheme results in describing the various singularity types by a pair (Si, Nj), where the first 
component reflects the behaviour of the Hubble parameter according to Theorem 13.2.14 
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while the second component symbolises one of the following three possible behaviours of 
the scale factor on approach to the finite-time singularity at t s namely, 

Ni a — > -collapse singularity or 

N 2 a — > a s 7^ -sudden singularity or 

N 3 a — > oo -hp singularity. 

We thus end up with nine possible singularity types given by the pairs (Si, Nj) where 
i,j = 1,2,3. We further note that in the (S2,N 3 ) category we can accommodate singu- 
larities having a Hubble rate which either diverges more rapidly or more slowly than the 
scale factor, as well as singularities with these two functions diverging almost or exactly 
at the same rate. Making such a further or more refined distinction between singularities 
that otherwise belong to the same (S, N) type means that we allow for the relative be- 
haviours of the Hubble rate and the scale factor to emerge as we approach the singularity. 
This refinement can be accomplished by introducing in our analysis the following notion 
of asymptotic strength. 

First recall that for the relative asymptotic behaviour of two functions /, g on approach 
to a point t*, we say that: 

1. f(t) is much smaller than g(t), f(t) « g(t), if and only if 
lim^ t . f(t)/g(t) = 

2. f(t) is similar to g(t), f(t) ~ g(t), if and only if ^ lim t _> t „ f(t)/g(t) < oo 

3. f(t) is asymptotic to g(t), f(t) <-> g(t), if and only if lim t _ t , f(t)/g(t) = 1, 

There are three specific (S, N) pairs that can be further analysed according to the 
possible relative asymptotic behaviours of H and a: 

We note that while (Si,N 3 ) and (S 2 , N 3 ) singularities can accommodate all possible 
relative behaviours between the Hubble rate and the scale factor, an (S3, N 2 ) singularity 
allows only for an asymptotic or similar relative behaviour between them. 
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We have so far managed to derive a new classification that consists of nine types of 
singularities described by the pairs (Si, Nj) where the first component reflects the behaviour 
of the Hubble rate (according to Theorem 13 .2 .11) . while the second one traces the behaviour 
of the scale factor. This new classification can also be further refined with the use of the 
notion of asymptotic strength. 

3.4 Bel-Robinson energy and singularities 

Even though we have at this point taken care of the two main functions that describe a 
singularity, we still have to face another aspect of the problem of classifying singularities. 
We have not yet taken into account the relative behaviour of the various matter components 
as we approach the time singularities. A clear picture of the asymptotic behaviour of the 
matter fields will provide us with the physical qualities of the singularity contributing 
significantly to the understanding of its nature. In order to describe the behaviour of the 
various matter components and include it in our geometric approach, we now exploit the 
use of the Bel-Robinson energy by combining the behaviour of the electric parts of the 
Bel- Robinson energy with that of the Hubble parameter in the following way: Suppose for 
example that S2 holds. Then from Eq. (IC.l.lOP we see that for a flat or closed universe \D\ 
necessarily becomes divergent whereas for an open universe it can be finite since the term 
k/a 2 may, in the case of some collapse singularities, exactly counter-balance the divergence 
of the H 2 term. These different cases can in turn be combined with a finite or an infinite 
\E\. Similarly, we can deduce the possible behaviours of \E\ and \D\ for a singular universe 
satisfying S\ or S3. Thus the necessary conditions for singularities based on the behaviour 
of the Hubble parameter given by Theorem 13.2.11 can now be rephrased in terms of the 
electric parts of Bel-Robinson energy, to read: 



Singularity type 



Relative asymptotic strength 



(Si,iV 3 ) 
(S 2 ,N 3 ) 
(S 3 ,N 2 ) 



H >> a, H « a, H <-> a, H ~ a 
H » a, H « a, H <-> a, H ~ a 
H <-> a, H ~ a 
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Theorem 3.4.1 Necessary conditions for null and timelike geodesically incomplete globally 
hyperbolic and regularly sliced isotropic universes are that at a finite time: 

Bi \E\ — > oo and \D\ — > oo 

B 2 \E\ < oo and \D\ — > oo 

B3 \E\ — > 00 and |D| < 00 

B 4 l^l < 00 and |D| < 00 

We can now move on to list all possible types of singularities that are formed in an isotropic 
universe during its cosmic evolution and enumerate the possible types that result from the 
different combinations of the three main functions in the problem, namely, the scale factor 
a, the Hubble expansion rate H and the Bel-Robinson energy B. These types will by 
necessity entail a possible blow up in the functions |J5|, \D\. 

If we suppose that a given model has a finite time singularity at t = t s , then the possible 
behaviours of the functions in the triplet (H,a, (\E\,\D\)) in accordance with Theorems 
I5XT1 13XT1 are as follows: 

51 H non-integrable on [t\, t] for every t > t\ 

5 2 H — > 00 at t s > t\ 

53 H otherwise pathological 

N 2 a — >• a s 7^ 

N 3 a — > 00 

B x \E\ — * 00, \D\ — > 00 

B 2 \E\ < 00, \D\ — > 00 
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B 3 \E\ -> oo, \D\ < oo 
B 4 \E\ < oo, |£>| < oo. 

The nature of any prescribed singularity is thus described completely by specifying the 
components in a triplet of the form 

(S l ,N„B l ), 

with the indices i,j,l taking their respective values as above. 

Note that there are a few types that cannot occur. For instance, we cannot have an 
(S 2 , N 2 , B 3 ) singularity because that would imply having a < oo (N 2 ) and H — > oo (S 2 ), 
while 3 ((a/a) 2 + k/a 2 ) 2 < oo (B 3 ), at t s which is impossible since \D\ 2 — > oo at t s (k 
arbitrary) . 

A complete list of the impossible singularities is model-dependent and is generally 
given by triplets (S i: Nj, B{) where the indices in the case of a flat or a closed universe take 
the values i — 1,2, j — 1, 2,3,/ = 3, 4, whereas in the case of an open universe they take 
the values % — 1, 2, j — 2, 3, I — 3, 4 (here by Si we denote for simplicity only the big bang 
case in the Si category). We thus see that some singularities which are impossible for a 
flat or a closed universe become possible for an open universe. Consider for example the 
triplet (S 2 , Ni, B 3 ) which means having H — > oo, a — > and 

|/J| 2 = 3((d/a) 2 + A;/a 2 ) 2 <oo 

at t s . This behaviour is valid only for some cases of an open universe. With these ex- 
ceptions, all other types of finite time singularities can in principle be formed during the 
evolution of isotropic matter-filled models in general relativity or other metric theories of 
gravity. 

This classification allows us to speak consistently of strong, mild or weak singularities. 
Note that the standard dust or radiation-filled big bang singularities fall under the strongest 
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singularity type, namely, the type (Si, Ni, Bi), since for example, in a flat universe filled 
with dust we have 

a(t) oc t 2/3 ^0, (Ni), (3.4.1) 
H oc r 1 -> oo, (3.4.2) 
\E\ 2 = 3/4# 4 ^oo, |Z?| 2 = 3H 4 -> oo, (Bi). (3.4.3) 

Our classification scheme is organized in such a way that the character of the singularities 
(i.e., the behaviour of the defining functions) becomes milder as the indices of S, N and 
B increase. Milder singularities in isotropic universes are thus expected to occur as one 
proceeds down the singularity list. 

It is the purpose of this classification to apply both to vacuum as well as matter 
dominated models. In fact, the Bel- Robinson energy takes care in a neat way the matter 
case. For instance, in fluid-filled models the various behaviours Bi, i = 1,2,3,4, of the 
Bel-Robinson energy density can be instead given by asymptotic conditions imposed on 
the density and pressure of the cosmological fluid: 

B 1 /i — ► oo and \/i + 3p\ — > oo 

Bi <=>■ /i — > oo and \fi + 3p\ < oo 

B 3 (j, < oo and \fi + 3p\ — > oo <^ ji < oo and \p\ — > oo 
B A <^ /i < oo and \ji + 3p| < oo <^ /i < oo and |p| < oo. 

We can translate these conditions to asymptotic behaviours in terms of a, H, depending 
on the value of k, for example, 

1. If k = 0, ji < oo =3- H 2 < oo, a arbitrary 

2. If k = 1, n < oo ^> H 2 < oo and a^O 

3. If k — -1, /j, < oo ^ H 2 - 1/a 2 < oo. 
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It is useful at this point to name five particular members of our classification of singu- 
larities that are of special interest. 

Definition 3.4.1 Let (Si, Nj, B k ) be a triplet describing the type of a singularity. 
We call the triplet: 

i) (S 1 ,N 1 ,B i ) ((S 2 , Ni, Bi)) a past-collapse (future-collapse) singularity. 

ii) (S 2 , N 3 , Bi) a big-rip singularity, 
ii) (S3, N 2 , -B3) a sudden singularity. 

iv) (S3, N 2 , Bi) a sudden-rip singularity. 

These types of singularities arise in some of the most popular cosmologies appearing in the 
recent literature (a detailed analysis of such cosmologies is given in Chapter 4). 

In order to gain a deeper insight of the dynamics on approach to the singularity 
we can make use of the notion of asymptotic strength defined in the previous Section. 
The use of asymptotic strength apart from offering us a clear view of the way that the 
three main functions, designating the nature of a singularity, compete asymptotically, and 
helping us decide about how strong one singularity type is compared to another, it also 
provides us with a refinement since it distinguishes singularities that fall in the same type. 
For example, as we have mentioned already all standard big-bang type singularities met 
in the literature belong to our past-collapse type (Si, N ± , Bi). Finding and comparing 
the asymptotic strengths of the singularities that may be present in dust and radiation 
filled universes we find that their precise nature is different since the radiation ones are 
characterised by the asymptotic strength a « H « (\E\ <-> \D\), whereas the dust ones 
have a << H « (\E\ ~ We shall give many examples of this later in Chapter 4. 

3.5 Closed trapped surfaces and Bel- Robinson energy 

In this section, we focus on the possible implications of the existence of a closed trapped 
surface in an isotropic universe and show how its existence can be traced in the behaviour of 
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the Bel-Robinson energy. We also illustrate through various examples how this association 
can be used to shed some light on certain aspects of the dynamical evolution of isotropic 
models. 

A closed trapped surface is a two dimensional spherical surface S in spacetime with 
the characteristic property that the expansion 9 of both the incoming and outgoing bundle 
of null geodesies orthogonal to S is negative as evaluated on S. That is, the light rays 
that are emitted perpendicularly from and towards S converge. This is feasible when the 
surface S encloses enough matter to attract all the incoming and outgoing rays of light. 
Since nothing can travel faster than light, the matter inside S becomes trapped within 
the boundary of this surface. If moreover the energy conditions hold (gravity remains 
attractive), the radius of S will decrease gradually until a singularity is formed pQ p. 2-3. 

In particular, sufficient conditions for the formation of a singularity are given by the 
well-known singularity theorems and include the existence of a closed trapped surface 
as well as appropriate energy and causality conditions imposed on the geometrical and 
topological structure of the spacetime (cf. [T|, p. 263-267): 

Theorem 3.5.1 Spacetime (V,g) cannot be null geodesically complete if: 

(1) R ab K a K h > for all null vectors K a 

(2) there is a non-compact Cauchy surface in V 

(3) there is a closed trapped surface in V. 

Theorem 3.5.2 Spacetime (V,g) is not timelike and null geodesically complete if: 

(1) R a t,K a K b > for every non-spacelike vector K a 

(2) every non-spacelike geodesic contains a point at which K[ a Rb] c d[ e Kf]K c K d ^ 0, where 

K a is the tangent vector to the geodesic 

(3) the chronology condition holds on V (i.e. there are no closed timelike curves) 

(4) there exists a closed trapped surface. 
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We see that the formation of a closed trapped surface is triggered by the presence of a 
strong gravitational field. A situation where this is possible is during the gravitational 
collapse of a star. In particular, when a star of mass twice that of the sun has burned 
out all the nuclear fuel which was supporting it against the influence of its own gravity 
and starts collapsing, a singularity becomes inevitable after a closed trapped surface is 
formed. A closed trapped surface is believed to have occurred even in the past of the 
universe itself. Looking back in time the universe enclosed enough matter to provoke a 
past trapped surface that led to a past-collapse singularity in the beginning of the universe. 

3.5.1 Existence 

Consider a past-directed bundle of radial null geodesies in an isotropic universe with metric 



where v is an affine parameter along these curves. The geodesic equation 

d 2 x s dx a dx b 

~d^ + Tab l^l^ = ^ (3 ' 
with x s = (t, r, 8, 0), for this type of geodesies has the following non-zero components: 



given by Eqs. (12X2) . (12X3) . Then 



2 d6 d(f) 
ds = — = — = 0, 

dv dv 



(3.5.1) 




0. 



(3.5.3) 




d 2 r dr dt 



0. 



(3.5.4) 



dv 2 r dv dv 



After substitution of the Christoffel symbols we obtain 




0. 



(3.5.5) 




d(t) dr dt 



a(t) dv dv 



0. 



(3.5.6) 
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Eq. (13.5.61) can be written in the form 



d 2 r/dv 2 
dr/dv 



dv 



and integrated out to give 



dr 



a(t) 



C 



(3.5.7) 



(3.5.8) 



dv a?{ty 

where C is an arbitrary constant and the sign depends on whether the null geodesies are 

incoming (+), or outgoing (— ). Inserting this in Eq. (13.5.51) (with C — 1), multiplying by 
dt/dv and integrating, we find 

dt 1 . 

— = ±— -, 3.5.9 

dv a(ty v ; 

from which we choose the minus sign since we are looking for the past directed family of 
null geodesies. The resulting tangent vector field is then 



k a 



dx a 
dv 



which gives 



dx a dx b 



9ab- 



a{t) ,± a 2 {t) , ° , ° 



s) +*«>(=) 



dv dv 

as desired. The divergence of the family of geodesies is calculated to be 



(3.5.10) 



(3.5.11) 



k a „ 



i.e., 



1 d_ 

'—g dx° 



k a „ 



-a(t) ± 



f(r) 



(3.5.12) 



(3.5.13) 



a\t) V v ' f(r) 

where f(r) = sin r,r, or sinhr for k = 1,0, or —1 respectively, [91] . A closed trapped 
surface is then formed if the divergence of both families of null geodesies becomes negative, 
that is if 

f(r) 



a(t) > 



f(r) 



(3.5.14) 
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for some values of r, t. 

An interesting fact for our purposes is that the above condition for the existence of 
a closed trapped surface has a decisive effect on the electric part \D\ of the Bel- Robinson 
energy. This is seen clearly if we use Eqs. f lCl.9p - flCl.10p to write \D\ in the form 



\m ^ 



a*(t) 



/(r)A W fa)) f 2 (r) 



(3.5.15) 



and then use the closed trapped surface condition (13.5.141) . In this case, the relation (I3.5.14p 
becomes equivalent to the inequality 

W > (3.5.16) 

from which it follows that collapse singularities (as predicted by the existence of a trapped 
surface) are characterised by a divergent Bel-Robinson energy. 

In the next Subsection we illustrate how some of the results found here help to deter- 
mine the evolution of various model universes. 



3.5.2 Examples 

Our first three examples illustrate the result that in the case when the CTS condition is 
combined with the null energy condition in spatially flat universes, an initial or future 
singularity is predicted to occur and this is reflected in the behaviour of the Bel-Robinson 
energy. In contrast, our last example 13.5.41 shows that even though the existence of a 
closed trapped surface is usually associated with the formation of a singularity this need 
not always be the case. This was recently shown by Ellis in [HI] who pointed out that a 
closed universe can escape the initial singularity provided that the non-spacelike energy 
condition does not hold. The null energy condition jointly with the condition for the 
existence of a closed trapped surface are not sufficient to imply a singularity in the case of 
a closed universe since Theorem 13.5.11 applies only in the case of an open or a flat universe. 
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Example 3.5.1 (Graduated inflation) 

Consider the flat graduated inflationary model developed in [92J, a universe filled with a 
fluid with an equation of state fi + p = 7/i 3 / 4 and 7 > (for more details on this family of 
models see Subsection 4.1.3). The null energy condition.!] is then automatically satisfied. 
The scale factor for this universe is given by 

fl W = ex P (-337^) ' ( 3 ' 5 - 17 ) 

(this is derived by setting 7 = 3/4 in Eq. (14.2.231) of Subsection 4.1.3), and the inequality 
Eq. (I3.5.14p for this case is equivalent to 

16 1 / 16 \ 1 , . 

= T7WTI? ex P ~~oTF? 57 > -• ( 3 - 5 - 18 ) 



3 3/2 7 2 f 2 V y 3 3/2 7 2 t J " r 

Hence we find that the condition for a CTS becomes valid for large values of r leading to 
the existence of past CTSs. Therefore the assumptions of the singularity theorem 13.5.11 are 
fulfilled making this model null geodesically incomplete in the past. 

Example 3.5.2 (Massless scalar field) 

A similar type of singularity is studied in [93J where we meet a flat model filled with a 
massless scalar field (see Subsection 4.1.1). If we consider a scalar field with <ft 2 > then 
the null energy condition is satisfied. The inequality (13.5. 14[) for this case reads 

a(t) = ^>K (3-5.19) 

which is again true for large values of r. Then using the singularity theorem 13.5.11 we 
conclude that this model is singular, and by the Bel-Robinson energy technique we arrive 
at a past-collapse singularity of type (Si, Ni, Bi). 

Example 3.5.3 (Sudden singularity) 



1 This condition leads to Eq. (|C.1.14[) which is equivalent to fj, + p > for the case of a fluid of density 
H and pressure p. 



3.5. CLOSED TRAPPED SURFACES AND BEL-ROBINSON ENERGY 



35 



A different situation is described by Barrow in [56] where he considers a model with scale 
factor 



<*(*) =[ T J (a. - !) + 1 - - T J , (3-5-20) 

with l<n<2,0<g<l and t s arbitrary but fixed from the beginning. This model 
experiences a sudden singularity (see Definition 13.4. II ii) and Section 4.4) at the future time 
t s with a(t) — > a^a(t) — > d s and d(t) — > — oo. For this model we find that the non-spacelike 
energy condition_| is valid because of a(t) — ► — oo, while for the case of a closed universe 
there exists a future CTS (the analogous inequality of (13.5. 14[) is a(t) < — \f'(r)/f(r)\) 
since the relation 

Q Tl 

<*(*) = (°« - !) + ^(*- ~ t ) n_1 < "I cotr l' ( 3 - 5 - 21 ) 
tl 

is satisfied for r = tt/2, a s < 1 and < £ < £ s . The singularity theorem 13.5.21 now applies 
and guarantees that this model is timelike and null geodesically incomplete in the future, 
while \E\ — > oo and \D\ < oo. 

Example 3.5.4 (De Sitter universe) 

The possibility of evading the singularity is feasible in a de Sitter universe for which a(t) = 
A cosh if t, k — 1, f(r) = sinr, so that the existence condition given by Eq. (13.5. 14[) is 
satisfied if and only if 

AHsinhHt > | cotr|, (3.5.22) 

which holds true for any t > and r = n /2 [91]. This implies the existence of more than 
one CTSs: every 2-sphere (t,r = const.) with t > and area greater than 

A s = An A 2 cosh 2 Ht sin 2 (r s ), (3.5.23) 



2 As we saw for a null vector this condition is equivalent to fx + p > 0. In order to find an analogous 
inequality for a timelike vector we multiply the Einstein equations by a unit timelikc vector k a : R a i > k a k b = 
{Tab - {l/2)Tg ab )k a k b = T ab k a k b + (1/2)T. Then R ab k a k b > <f> T ab k a k b > -1/2T. For a perfect fluid 
this translates to \x + 3p = —6d/a > 0. Therefore the non-spacelike energy condition holds iff /i + 3p > 
and /i + p > 0. 
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where cot r s = AH sinh Ht, is a CTS. Nevertheless the spacetime is complete. The timelike 
energy condition is violated, (/i + 3p) = —d(t)/a(t) < 0, and so Theorem 13.5.21 does not 
apply. However, the null energy condition (1C. 1.121) holds true but Theorem 13.5.11 does not 
apply because the spatial slices are now compact. Then we find that since 



\E\ = V3 



a(t) 



a(t) 



V3H 2 



(3.5.24) 



and 



\D\ = V3 



a 2 (t) + 1 



a 2 (t) a 2 (t) 



V3 



H 2 tanh 2 Ht + 



(3.5.25) 



A 2 cosh 2 Ht ' 

they always remain finite indicating the complete character of the spacetime in this case 
(note that the assumptions of the completeness theorem of [82] which we discussed in the 
end of Section 2.2 are satisfied). 



3.6 Discussion 

We started off in this chapter by giving a classification of the singularities that can occur 
in an isotropic universe based on the behaviour of the Hubble parameter. However, such 
a classification cannot lead to complete results and we have found it necessary to extend 
this scheme by including succeedingly the behaviours of the scale factor and of the Bel- 
Robinson energy. In this case, we have found that the resulting behaviours of the three 
functions, H, a and B taken together exhaust the types of singularities that are possible 
to form during the evolution of an isotropic universe. 

The resulting classification is described by triplets of the form (Si, Nj, B{) where the 
S category monitors the asymptotic behaviour of the expansion rate, closely related to the 
extrinsic curvature of the spatial slices, the N category that of the scale factor, describing 
in a sense what the whole of space eventually does, while the B category describes how 
the various matter fields contribute to the evolution of the geometry on approach to the 
singularity. We know (cf. [72], [94]) that all these quantities need to be uniformly bounded 
to produce geodesically complete universes. Otherwise, the whole situation can be very 
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complicated and we have exploited what can happen in such a case in the relatively simple 
geometry considered here. 

Our scheme not only covers all the recently discovered types of singularities but it also 
predicts many possible new ones. For example, in the case of a flat isotropic universe the 
classification of [52] provides us with four main types of singularities. These four types 
can be identified with four particular (Si, Nj, Bk) triplets of our scheme: Namely the 'big 
rip' type of [52] characterised by a — > oo, \i — > oo and \p\ — > oo at t s is an (S2,N 3 ,Bi) 
singularity; their 'sudden' singularity described by a — > a s < oo, p, — > fi s < oo and \p\ — > oo 
at t s is an (S 3 , N 2 , B 3 ); further, their type III singularity, namely a — > a s < oo, p, — > oo 
and \p\ — ► oo at t s , is clearly an (5*2, N 2 , Bi) type singularity, while type IV singularities of 
[52] with a — > a s < oo, /i — > and \p\ — > at t a , all belong to the (5*3, A^ 2 , -B4) category. 

Having all possible singularity types in the unified form expounded in this chapter has 
the added advantage that we can consistently compare the various different types as we 
asymptotically approach the time singularity. In this case we saw in detail that the relative 
strength of the functions describing the singularity type becomes an important factor in 
order to distinguish between all possible behaviours. 

It is natural to consider the extension of our classification scheme in the context of 
more general classes of cosmological spacetimes. We believe that an analysis of this more 
complicated case is still feasible using the techniques of the present work. In particular, 
the first step towards that direction would be to examine the family of homogeneous but 
anisotropic universes. The nature of a singularity arising in such a universe would be given 
by a triplet analogous to that of the isotropic case but now the number of possible types of 
singularities will increase. This will result from the existence of more than one scale factors 
describing the model and also from the fact that the magnetic parts of the Bel-Robinson 
energy will not be zero as they were in the isotropic case studied here. 

A further step in that direction would be to examine inhomogeneous spacetimes where 
the space gradient of the lapse will be a function not only of the time but also of the space 
variables and hence it will not vanish automatically. The nature of a singularity would 
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then be given by a quadruplet with the new component describing the possible behaviours 
of the space gradient of the lapse. These behaviours are predicted from the contrapositive 
of the completeness t heorem 12 . 2 . 1 1 and they include both a diverging | VN\ g in a finite time 
(giving rise to a possible blow-up lapse singularity) as well as a |ViV| s which will be finite 
only for a finite interval of proper time corresponding to a sudden lapse singularity. The 
analysis of all these cases lies beyond the work in the present thesis. 



Chapter 4 

Singularities in isotropic cosmologies 



4.1 Introduction 

In the previous chapter we derived a complete classification of spacetime singularities that 
are in principle possible to arise in isotropic cosmologies. The purpose of this chapter is 
to look for physical circumstances that give rise to four particularly important members 
of our classification given by Definition I3.4.H ideally we would be content with prototypes 
of cosmologies with every possible singularity type that is predicted by our classification 
scheme. 

The analysis of the various cosmologies appearing in this chapter is organised as fol- 
lows. In Section 4.2 we study collapse singularities that arise in recently studied cosmologies 
containing a massless scalar field (Subsection 4.2.1), a scalar filed with an exponential po- 
tential (Subsection 4.2.2), or a graduated inflationary potential (Subsection 4.2.3). Section 
4.3 provides an analysis of big-rip singularities which are met in models containing dark 
energy in the form of phantom fields (Subsection 4.3.1), or Chaplygin gases (Subsection 
4.3.2). This type of singularity offers a new possibility of a violent end for the universe 
towards which every bound existing structure gets 'ripped apart'. In Section 4.4, we ex- 
plore the character of sudden singularities that may arise in cosmologies filled with a fluid 
having either unconnected pressure and density (Subsection 4.4.1), or alternatively, obey- 
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ing an equation of state (Subsection 4.4.2). In Subsection 4.4.3, we analyse models that 
exhibit milder than sudden singularities. Finally, in Section 4.5, we discuss sudden-rip type 
singularities. 

In each one of these sections we prove theorems giving necessary and sufficient condi- 
tions for the occurrence of the various types of singularity that are analysed. We expect the 
proofs of all these theorems to be quite straightforward for we have now already identified 
the type of singularity that we are looking for in accordance to our classification. Proving 
such theorems without this knowledge would have been a problem of quite a different order. 

The usefulness of these theorems lies in that they answer the question of whether or 
not the behaviours met in known cosmological models described by exact solutions (which 
as a rule have a flat spatial metric (k = 0)) continue to be valid in universes having nonzero 
values of k, or described by solutions which are more general than being exact in the sense 
that some or all of the arbitrary constants present in a solution remain arbitrary. We 
will see that whenever the behaviour met in flat models passes over to curved ones either 
the curvature term in the Friedman equation turns out to be subdominant compared to 
the density term, or, in any case, cannot significantly alter the behaviour of the Hubble 
parameter H. 

It is worth mentioning that all the different example cases of cosmological models 
analysed in this Chapter have been a primary source of motivation for finding more general 
theoretical criteria which determine their evolution and mark their basic characteristics. 

4.2 Collapse singularities 

Following the Definition 13.4.11 given in the previous chapter, by past or future collapse 
singularities we mean those that arise in solutions with a scale factor tending to zero at an 
initial time (big-bang time), or within a finite future time (big-crunch time) respectively, 
while the Hubble expansion rate and the electric parts of the Bel- Robinson energy diverge. 
Originally, collapse singularities were discovered in simple FRW universes filled with dust 
or radiation (cf. [77j, pp. 96-101). Here, we will instead focus on recently analysed models 
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that involve the presence of a scalar field (cf. [93], [95], [92], [96J). 

Scalar fields play a substantial role in cosmology. In the framework of inflation, scalar 
fields having a self-interaction potential are used to model the matter content of the early 
universe and their dynamical evolution monitors the mechanism of inflation (cf., e.g., [9~Tj . 
pp. 261-317). However, the role of scalar fields is not limited to early universe scenarios. 
They are also crucial in the present phase of evolution of the universe since scalar fields 
actualise a possible dark-energy component that is believed to drive the presently observed 
accelerating expansion of the universe (cf. [EE]). Furthermore, in grand unified theories 
scalar fields are met in the spontaneous symmetry breaking (cf. [97], pp. 195-255, |98j ) 
and other fundamental mechanisms proposed within these theories (cf. |97j). 

In a homogeneous universe, scalar fields are functions of the proper time only. The 
energy-momentum tensor of a scalar field is equivalent to that of a perfect fluid with density 
/x^ and pressure p^ given by (cf. Appendix D), 

The relation between the pressure and the energy density of the scalar field is therefore 

P<t> = H~ 2V > 

which is not an equation of state in the usual sense [99] . 

However, in the special case of a massless scalar field, V(4>) = 0, we have an equation of 
state of the form p$ = il<$>, that is it behaves as stiff matter. If on the other hand, the energy 
of the scalar field is purely potential, or dominates over the kinetic energy, then the scalar 
field is approximately described by the equation of state = —li^, and behaves therefore 
like a cosmological constant. Within the framework of inflationary cosmology for example, 
the energy of the scalar field becomes mainly potential as soon as the field enters the 
slow-roll phase. This means that <j) 2 <C V(<fi), or equivalently /x^ + 3p<^ < 0, and hence the 
dominant energy condition is violated in this case. Therefore the singularity theorem 13.5. II 
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does not apply and the initial singularity is not an inevitable consequence of inflationary 
cosmologies. Nevertheless, there exist interesting, non-trivial scalar field cosmologies that 
cannot avoid the initial singularity in this way and in the next subsections we shall focus 
on the analysis of such cosmologies. 



4.2.1 Massless scalar field 

An example of a singular scalar field cosmology that possesses a past-collapse singularity 
can be found in |93j. This model is a flat FRW universe with a massless scalar field, and 
the Friedman and continuity equations read 

3H 2 = l<p 2 , + 3#0 = O, (4.2.2) 



from which we readily find: 



Indeed, this universe emerges from a past-collapse singularity at t — 0, since a oc t 1 / 3 - 0, 
H -> oo, while \E\ 2 = 4 /3 -> oo and |£>| 2 = 3H 4 -> oo. 

In [93] Foster also found a condition that should be satisfied so that solution (14.2.31) . 
approximates that with an arbitrary potential V(4>) in the limit t — > 0. This condition is 
derived by simply substituting the solution given by Eq. f)4.2.3[) in the Friedmann equation 
that is satisfied by a scalar field with potential V(<p), i.e. 3H 2 = V(<f>) +4> 2 /2, and equating 
the resulting terms. One then finds that 

lim t v(±y|ln^=0, (4.2.4) 

or 

lim e-^V^) = 0. (4.2.5) 

cj>— >±oo 

Hence the solution given by Eq. ( 14.2.3ft approximates the solution for a scalar field with 
potential V((f>) when the potential decreases slower than e - ^'^'. The following theorem 
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implies that this type of singularity may also arise in more general, curved universes [69] . 



Theorem 4.2.1 A necessary and sufficient condition for the occurrence of a past-collapse 
singularity of type (S%,Ni,Bi) at time t x in an isotropic universe with a massless scalar 
field is that ^ oo att\. 

Proof. We first prove necessity. Integrating the continuity equation + 3H<fi = 0, we 
find that oc a~ 3 . Since —* oo, we have that a — > 0. Because 



H 2 = ^- — ^oo, as (4.2.6) 
3 er 

H becomes unbounded at t\. In addition, since 

1 1 3 12 v ; 



and 



|£| 2 = ^(M + 3p) 2 = ^-oo, (4.2.8) 



as t — *■ ti, both \D\ and \E\ diverge there. Conversely, assuming a past-collapse singularity 
of type (Si, Ni, Bi) at ti, we have from condition Bi (<$ \i — ■> oo, |/i + 3p| — > oo) that 
/x — > oo and so 2 — > oo as t — > ti. ■ 

It follows that the asymptotic strength of the singularity in this model (described by 
both the exact solution and the above theorem) is given by: a « H « (\E\ ~ \D\). 



4.2.2 Multiple exponential potential 

Since there is not as yet any unique well-established fundamental theory dictating and 
predicting an exact form of the potential V(0), there are many possible candidate potentials 
currently studied depending on the underlying theory they intend to describe. For example, 
in [95] the authors consider a scalar field with the 'multiple' exponential potential of the 
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V(<f ) ) = W -V smh(^3/2<j ) ), 

where Wo and Vq are arbitrary constants The scalar field inhabits a flat FRW universe 
containing also a separately conserved pressureless fluid. In this context, the field equations 
become 

3H 2 = (i+V(<f>)+v), (4.2.9) 

2H + 3H 2 = - (t^-V(<f>)\ , (4.2.10) 

4> + 3H<j) + V'((j)) = 0, (4.2.11) 
fi + 3Hfx = 0. (4.2.12) 

Following [95], we split the scale factor in a product form by setting a 3 = xy, while 
the scalar field is written as = a/2/3 \og(y/x) with a; and y positive. Setting 

x = C[expxicosx 2 + exp(-xi) C0SX2], (4.2.13) 
y = C[expxisinx 2 + exp(-xi) sin(-x 2 )] (4.2.14) 

(here, Xi — w i(t — to)iX2 = w 2 (t — T ), C > 0, t is an arbitrary constant, T is the 'initial' 
time and W\, w 2 positive parameters such that w\ — w\ = 3/4W ; 2wiW 2 = 3/4Vq), we 
find that for large and positive values of the time parameter t the scale factor becomes 

a = C 2/3 exp(2/3xi) cos 1/3 \2 sin 1/3 X2, (4.2.15) 

which is obviously divergent. 

At the finite time value T s = 7r/(2w 2 ) + T , we have that a — > 0, H — > u^2(cotx 2 — 
tanx 2 )/3 — > 00, and consequently we find that 



(4.2.16) 
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and also 

\E\ -> |2w 2 (2(cotx2 - tanx 2 ) - ^(cot 2 x 2 + tan 2 x 2 ))| -> oo. (4.2.17) 

9 

This model therefore serves as a nice example that illustrates the structure of the future- 
collapse type singularity predicted earlier in our classification. 



4.2.3 Graduated inflation 

In the original inflationary scenario the equation of state p = —fi that describes the scalar 
field during the slow-roll period is eventually responsible for an accelerating exponential 
expansion. However, milder forms of solutions, for example a power-law expansion a (xt p 
with p > 1, can still suffice for inflation. Different forms of solutions in which the expansion 
follows a power-law with an arbitrary in general exponent lead to the so-called intermediate 
or graduated-inflationary universes (cf. [92] and references therein). These models were 
first studied by Barrow who introduced an equation of state of the form 

H+p = jfi\ (4.2.18) 



where A, 7 are constants and 7 7^ 0, in a flat model. The set of solutions that Barrow found 
(which naturally depended on the A and 7 parameters) included not only the traditional 
power-law and exponential inflationary solutions but also new ones. In these solutions the 
scale factor is an exponential of an arbitrary power of the time and in some cases describes 
a universe that begins with a past-collapse singularity. In the original references, however, 
no complete discussion of the structure of singularities predicted in these models is included 
and so we give below such an analysis. 

To see how a past-collapse singularity develops in the history of these cosmologies we 
substitute the equation of state (14.2.181) in the continuity equation and integrate to find 

(37(A - l))V(i-A) 
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where A 7^ 1 and a$ is a constant. Using then the Friedmann equation 3H 2 = fi, we obtain 
the scale factor in the following explicit form [92J: 

/ fl \ oA/(l-2A) l/(l-2A)/x _ 1 /o\2(l-A)/(l-2A) 

In ( - ) = 1 f- i/Z) t 2(l-A)/(l-2A) (4 2 2Q) 

\ a 0/ A-l 

with A ^ 1, 1/2. Setting 

oA/(l-2A) l/(l-2A)/x _ 1 / 2 ^2(1-A)/(1-2A) 

A = 1 | A , (4.2.21) 

A — 1 

/ = 2(1-A)/(1-2A), (4.2.22) 

we can write the scale factor in the form 

a(t) oc exp(A^). (4.2.23) 

Considering the case that 1/2 < A < 1 and 7 > 0, we have that A < and / < 0. This 
means that a past-collapse singularity develops as t — > 0, since then a -> 0, = A/t-^- 1 -> 
00 and consequently 

|£)| 2 = 3# 4 -> 00, (4.2.24) 
and also \i — > 00, which in turn implies that 

|£| 2 = -L(-2/i + 3 7 /i A ) 2 -> 00. (4.2.25) 

The asymptotic strength of this singularity is therefore: a << (\E\ <-> |D|) << i7. Actually 
this type of behaviour is found to occur more generally due to the following result. 

Theorem 4.2.2 A necessary and sufficient condition for a past-collapse singularity 

(Si, N\, B{) to occur at time t\ in an open or flat universe filled with a fluid with equation 

of state p + fj, = 7/i A , 7 > and A < 1, is that fi — > 00 as t — >• t\. 



Proof. The continuity equation Eq. ( 12.1. 6ft gives directly the form of the scale factor, 
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a = a exp \—^—-^), (4.2.26) 



so that a — > as t — > t\. Since 



2 n k n 2 / -2/i A+l 



|£| 2 = l(-2^ + 3 7 //) 2 , 



3 ' 



we see that as t — > t\ H, \E\ and \D\ diverge provided that k = or k = —1. The converse 
is straightforward. ■ 

As it was pointed out in [92] the matter content in a graduated inflationary universe 
can be reinterpreted as a scalar field with the following potential: 

f 1 \ (/-!)// \ 2 / f _ i \ (/-2)// 



(4.2.27) 
where 

U(l-f)\ 1/2 

J 1 (4.2.28) 



/ 

and where / is given by Eq. ( 14.2.22ft . It is worth noting here that it is only recently that 
this model has regained interest since, as it was shown in [96], a specific choice for the / 
parameter offers a good fit to the WMAP observations [55] within the framework of the 
slow-roll approximation. 

In order to illustrate this result, let us briefly review the equations of the slow-roll 
approximation. For large the potential ( 14.2.27ft becomes [96] 



V(<f>) = ^^(2^)^, (4-2.29) 
where = {2Aj3t^) 1 ^ 2 and (3 = 4(/ _1 — 1). In the slow-roll approximation it is assumed 
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that cf) 2 <C V and <fi <C V restricting the field equations to the form (cf. |100j . [101] ) 



Inflation then occurs if 



3H 2 ~ V (4.2.30) 

3H<j) ~ -V. (4.2.31) 

- = # + # 2 >0^-^r<l. (4.2.32) 
a H 2 



Differentiating Eq. (14.2.301) with respect to t and using Eq. (14.2.311) we find that 

(V) 2 

H= -m* (4 - 2 - 33) 



so that 

while condition (14.2.321) holds if and only if the first slow-roll parameter e, defined by [100J, 
[TOT] 



1 fV' x 2 



*=2\v) ' (42 - 35) 
satisfies the condition 

e < 1. (4.2.36) 

If we further differentiate Eq. (I4.2.3ip and use Eq. (14.2.361) together with </> <C V, we 
conclude that the second slow-roll parameter r], defined by [100] . [101] 

V" 

V = -y, (4.2.37) 

satisfies the condition 

\v\ < I- (4.2.38) 

However, as it is stated in [102] , it is not necessary to demand that the derivative of an 
approximation equation is itself a valid approximation equation and so this last condi- 
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tion can be regarded as independent from the first one given by Eq. (14.2.361) . fn the 
intermediate-inflationary model of [96J these parameters are shown to be e = (3 2 /2<j) 2 and 
rj = (3(1 + (3/2)/<j) 2 . As it is further shown in [96], the resulting values of two very impor- 
tant spectral indices that are closely related to these slow-roll parameters come in good 
agreement with their corresponding values measured by observations, provided that (3 = 2, 
i.e, / = 2/3 and hence a(t) = exp(At 2 ^ 3 ). A graduated-inflationary model with these 
parameters and scale factor offers therefore an appropriate theoretical framework within 
which the data from the WMAP observations is explained. 

4.3 Big-rip singularities 

Let us now move on to the discussion of a less standard type of model, nevertheless one 
which is singular in a way predicted by our scheme. The notion of dark matter was 
introduced in cosmology in order to explain a series of observations in the area of galactic 
dynamics (cf. [103] pp. 417-424). These observations showed that the velocities of galaxies 
in clusters are so high and the luminous masses enclosed within the galaxies so low that 
these structures ought to have been ripped apart. A similar situation holds for the case 
of spiral galaxies. In particular, it was found in [104] that at a distance r outside the 
central luminous part of such a galaxy the rotational velocity v(r) does not follow the 
anticipated third law of Kepler v(r) oc r -1 / 2 but is instead measured to be constant. This 
leads to the counter-intuitive fact that the mass at the outskirts of such a galaxy grows 
proportionally to distance (cf. [103] pp. 46-47). Such a pattern of matter distribution 
could be attributed to the existence of a halo surrounding these galaxies and consisting of 
low-luminosity matter, thereby called dark matter, that contributes the missing amount of 
matter needed to fit the observed profile of these galaxies (see [103] pp. 45-53). 

Dark matter is also considered in another context. Observations strongly indicate 
that the value of Q, the ratio of the present density of the universe to the critical density 
fx c = 3Hq/StiG (H being the present value of H), is unity. However, the luminous mass 
of the universe is currently measured to have a value of flium only around 0.04 [105] . On 



50 



CHAPTER 4. SINGULARITIES IN ISOTROPIC COSMOLOGIES 



the other hand, measurements of the amount of dark matter that is believed to exist in 
clustered areas of the universe give flom — 0.26 [105J. This missing part, nearly 70% of 
the critical density, is contributed by an unclustered component the dark energy which 
has a ratio of pressure to density w — pj\x < —0.78 [105J. Taking this into account we 
see that dark matter and dark energy may play an important role in the understanding of 
the composition of the universe. Even though the nature of dark energy remains largely a 
mystery it is an issue worth examining as it offers us the opportunity to discover new and 
exciting possible fates of model universes with such forms of material content. 

The simplest source of dark energy can be described by simply postulating the existence 
of a cosmo logical constant A which satisfies a perfect fluid equation of state, p = wfi with 
w = — 1 [106J. However, the consideration of a cosmological constant A in this context leads 
to the so-called cosmological constant problem - the expected value of the cosmological 
constant from quantum field theories is some 120 orders of magnitude larger than the 
observed one [3]. Another candidate for dark energy is the quintessence, a scalar field 
which satisfies the equation of state for — 1 < w < — 1/3 (cf. |107j - [TT3] ). Because recent 
observations indicate that the universe is not only expanding but is also accelerating a lot 
of attention has been drawn to use such matter components (or even more exotic ones, cf. 

in an effort to explain the observed acceleration [54 J [55]. Note that in order to 
obtain an accelerating FRW universe filled with a perfect fluid equation of state p = wfi, 
it is necessary to have a value of the w parameter that is less than —1/3 as this clearly 
follows from the field equations (cf. I4.3.3P below). 

Cosmological observations do not exclude the fact that the w parameter be less than —1 
[3]. An entity satisfying such an exotic equation of state is dubbed phantom dark energy; 
this is a becoming name for the description of something that is physically unorthodox but 
which nevertheless seems to provide an explanation to the observations [53] . To begin with, 
phantom dark energy immediately violates the dominant energy condition since p + \i — 
fi(w + 1) < 0. This leaves room for the existence of wormholes and times machines (see (3] 
for a discussion). Another interesting albeit unphysical property of phantom dark energy 
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is that it exhibits a sound speed that it is both imaginary and superluminal [2]. 

Despite the strange properties of phantom-dark energy it is still an interesting type of 
matter field to study in cosmology (see Refs. [2]-[lZ], [53], for phantom filled cosmological 
models) since its presence influences the future of the universe through the development 
of a big-rip singularity as is shown in the next Subsection. In what follows, we will use the 
term phantom cosmology to refer to any model that includes a phantom- matter component. 



4.3.1 Phantom dark energy 

Consider the case of a flat FRW universe with metric ds 2 = —dt 2 + a 2 (t)dx 2 , filled with a 
fluid with equation of state p = wp. The Einstein field equations for such a universe are 



p = 3H 2 , 

ft + 3Hp(l + w) = 0, 

- = --//(l + 3i«). 

a o 



(4.3.1) 
(4.3.2) 
(4.3.3) 



Integrating the continuity equation we have 



H = Ca~ 3{w+1 \ 



(4.3.4) 



where C is the integration constant. For ordinary matter (w > 0), or for dark matter 
with w > — 1, the energy density decreases as the universe expands. In contrast, phantom 
dark energy has a density that increases as the expansion proceeds in an FRW universe 
(irrespectively of its spatial geometry). This implies that it will at first overcome and 
finally dominate over all components of ordinary matter. 

Solving the Friedmann equation (14.3. ip for a, we find [51] 



QQ 3 ( i + ,)/ 2+ 3(l+^)^ (t _ / 







3(l+u.) 



(4.3.5) 



52 



CHAPTER 4. SINGULARITIES IN ISOTROPIC COSMOLOGIES 



where A = C/3 and ao is an integration constant. The behaviour that would have been 
expected had we considered only ordinary matter is that a — > oo as t — > oo, and the 
Hubble distance H^ 1 to grow more rapidly than the scale factor (as it follows from Eq. 
(14.3.51) and Eq. (I4.3.7P below), which means that more and more galaxies would become 
visible asymptotically [51]. In contrast to this situation, in a phantom dark energy model 
a big-rip singularity arises within the finite value of time given by [51] 



Uip -t + - — ■ 3(H _ 1)/2 , (4-3.6) 



3VA(\w\ - l)a, 







and called the big-rip time. Note that t r i P depends on the initial size of the universe, ao, 
and the exact value of the w parameter in such a way that the bigger the values of ao and 
w, the briefer the existence of the phantom universe. Indeed, from Eq. (14.3.51) we see that 
the scale factor becomes divergent as t — > t rip . Also the Hubble rate and hence the norms 
the electric parts of the Bel-Robinson energy become diverge there since we can obtain the 
following forms: 



H = A 



a ^^ + 3{1+ : )VJ (t-t ) \ (4.3.7) 



\Ef = ^ 4 (l + 3w) 2 , (4.3.8) 
|£>| 2 = 3H 4 . (4.3.9) 



Eq. ( 14.3.7!) implies in particular that the Hubble distance H~ l now tends to zero as t — > t r i P . 
Combining this with the infinite growth of the scale factor as t — > t r i P , we see that this 
type of accelerated expansion forces galaxies to disappear beyond the cosmological horizon 
rendering the end of the universe dark and rather violent. The latter attribute follows from 
the fact that it can be shown (see [3]) that as the big-rip time approaches gravitationally 
bound structures such as galaxies, planetary systems or even atoms, become unbound due 
to the enormous repulsive pressure. 



The following theorem [69], [70] characterises the future singularity in phantom cos- 
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mologies irrespective of the value of the curvature k and implies that the singularities in 
such models can be milder than the standard all-encompassing future-collapse singulari- 
ties and have necessarily diverging pressure. This theorem generalises and unifies previous 
results (cf. [3], [IE], [H], [20], [5T], [52], P], [H], [33). 

Theorem 4.3.1 Necessary and sufficient conditions for a big-rip singularity (S2, N 3 , Bi) 
to occur at the finite future time t s in a universe filled with a fluid with equation of state 
p = wji, are that w < — 1 and \p\ — > 00 at t s . 

Proof. Substituting the equation of state p = wji in the continuity equation fi + 3if (/i + 
p) = 0, we have 

fi oc a- 3(w+1 \ (4.3.10) 
and so if w < — 1 and p blows up at t s , a also blows up at t s . Since 

H 2 = ^-^, \D\ 2 = ^-, \E\" = ^(l + 3wf, (4.3.11) 
3 a z 3 12 

we conclude that at t s , H, a, \D\ and \E\ are divergent. Conversely, assuming a big-rip 
singularity (S 2 , N 3 , Bi) at t s in an FRW universe with the equation of state p = wfi, we 
have from the (Bi) hypothesis that fi — > 00 at t s and so p also blows up at t s . Since a is 
divergent as well, we see from (14.3. lOj) that w < — 1. ■ 

It is interesting to note that the relative strength of this big-rip singularity is strongly 
influenced by the choice of the w parameter. In particular, there are three possible be- 
haviours depending on the given ranges of the w parameter: 

• if -4/3<w<-l, then H « (\E\ ~ \D\) « a, 

• if -5/3 < w < -4/3, we have H « a « (\E\ ~ \D\), whereas 

• if w < -5/3, we have a « H « (\E\ ~ 

We note that for all three ranges of w, \E\ and \D\ diverge similarly while H remains 
always much smaller compared to these two quantities. On the other hand, the scale factor 
diverges less strongly relatively to H and \E\, \D\ as w decreases below —1. 
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4.3.2 Chaplygin gas 

As we saw in the previous Section a phantom equation of state drives the universe towards 
a big-rip singularity. In this section we will see that a big-rip singularity can still arise in 
models with forms of dark energy different from phantoms such as the so-called generalised 
Chaplygin gases; these fields play a central role in dark energy cosmology and specifically 
in the unification of dark energy and dark matter (see Refs. |114j - [T3~8] ). 

For example, consider a flat FRW universe containing a dark-energy component which 
at late times satisfies a general equation of state of the form [32] 

p = -H i -Afi (4.3.12) 



where A and a are real parameters. Integrating the continuity equation 



fi d - 3HA^ = 0, (4.3.13) 
we find that the energy density fid of the dark component behaves as follows, 



a 



Hd = Hdo[l + 3A(1 - a) In — ] , (4.3.14) 



where A = Afi^ 1 > and /i^o is an integration constant. Note that 

l 



l-ct 



dto _ZA ND f 1 + 3A{1 _ a)ln f L] (4 .3. 15) 



da a \ ao 

so that the parameter A monitors the behaviour of the dark energy as the universe expands: 
if A > the density of the dark energy increases with expansion if A = 0, the density is 
constant, whereas if A < 0, the density decreases [48J. Recall that for the occurrence of 
a big-rip singularity the density needs to grow with expansion and hence for the study of 
such a singularity we may focus here in the case A > 0. 

Substituting in the Friedmann equation f)2.1.4p we find the scale factor in the following 
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implicit form 



l-2a 1— 2a 



n \ 2(l-a) f _ n \ 2(l-a) Q 

1 + 3,4(1 -a) hi — - l + 3A(l-a)ln— = -Ail - 2a)C(t - t d ). 

(4.3.16) 

Here B = a/ /ido/3, td is the time at which the dark energy becomes dominant and ad is 
the corresponding value of the scale factor ad = a(t d ). The Hubble expansion rate is then 
given by 

l-2a 

a A2(i- Q) 3 



// B\ ( 1 + 3A(1 - a) In ^ J + -A(l - 2a)C(t - t d ) ) (1.3.17) 

We now show that the future behaviour of this model is quite sensitive to the values of the 
parameter a. Although the case with a G (1/2, 1) leads to the development of a big-rip 
singularity in finite time, when a < 1/2 we are led instead a singularity-free universe. 
Indeed, taking a G (1/2, 1) we allow for the dark energy density, the scale factor and the 
Hubble rate H all to increase with time until they diverge at the finite time tf where 



(hi 



2(l-a) 



''°' j+ M fa- 1)c r M "-° |h w • < 4 - 3i8) 

In addition, the Bel-Robinson energy blows up there since the electric parts take the forms: 

\E\ 2 = l( 2/U + 3A/i a ) 2 , (4.3.19) 
\D\ 2 = 3H A . (4.3.20) 



If instead one assumes that a < 1/2, we see from Eq. (14.3. 17p that H is always finite 
(so that the Bel-Robinson energy is also finite) and therefore the model is geodesically 
complete. (The case with a > 1 leads to an (S 2 , N 2 , Bi) singularity since as it follows 
from (14.3. 14p fid and consequently pd, H, \E\ and \D\ all diverge at the finite value a s = 

aoe l/(3A(o-l))_) 

Summing up the results of this Section, the consideration of these unconventional dark 
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energy components in the framework of FRW cosmology leads to the possibility of a big-rip 
singularity. Even though Theorem 13.2.11 predicts such a singular behaviour in condition 
5*2, it is only very recently that the study of physically motivated models described here 
has offered sufficient physical grounds for the appearance of such a structure distinct from 
the standard future-collapse types considered so far (cf. [3], [H], [19], [50], [51], [52], [53] . 
[H], [38]). 



4.4 Sudden singularities 

In the previous Section, we saw that the violation of the dominant energy condition is 
sufficient for the occurrence of future finite-time singularities during the expansion of the 
universe. But is this also necessary? Can a finite-time singularity develop during the 
expansion of the universe without violating the dominant energy condition? If yes, then 
what are the qualitative differences between the resulting types of singularities? 

Barrow initiated a study of these questions in [56] . He found a new type of singularity 
that can arise in a universe where matter does not violate the dominant energy condition. 
Actually in his example matter satisfies even the strong-energy conditions fi > and 
+ 3p > 0. This fact implies that the violation of the dominant-energy condition is 
sufficient but not necessary for the occurrence of a finite-time singularity during expansion. 
Moreover, the nature of the resulting singularity is substantially different from the big-rip 
singularities discussed earlier. In particular, the finite-time singularity of Barrow is a 
sudden singularity according to our Definition 13.4.11 occurring at an arbitrary but fixed 
value of finite time t s at which the scale factor, the Hubble rate and consequently \D\, 
all remain finite but \E\, or equivalently the acceleration, diverges. This becomes possible 
since by the Raychaudhuri equation at t s we have 

a p He k 

- = -- - -, 4.4.1 

a 2 2 2a;?' v ; 

where a s = a(t s ), H s = H(t s ) > 0, which implies that a singularity with a/a — > — oo can 
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occur provided the pressure blows up at t s . 

In the following Subsections we analyse how a sudden singularity arises in universes 
filled with a fluid having fi and p unconnected (Subsection 4.4.1), or satisfying an equation 
of state (Subsection 4.4.2). In Subsection 4.4.3 we explore model universes that exhibit 
singularities that are milder than the sudden ones. 

4.4.1 Fluids with p and fi unconnected 

A special exact solution acquiring all the characteristics needed for the development of a 
sudden singularity is given in [56] and reads 

a(f) = l+Q^V.-l) + (l-£r, (4- 4 - 2 ) 

with l<n<2, 0<g<l. For the more detailed analysis of the properties of this solution 
given below we follow [67], [68] ■ This solution exists smoothly only on the interval (0,t s ), 
with a s and H s and therefore \D\ is finite at the right end. By setting A = q(a s — l)/tf, B = 
njt n sl we find that 

d(t) = At"" 1 + B(t s - 1)"- 1 , (4.4.3) 

which means that unless q — 1, a blows up as t — > 0, making H continuous only on (0, t s ). 
Also a(0) is finite and we can extend H and define it to be finite also at 0, H(0) = Ho, 
so that H is defined on [0,£ s ]. However, since lim t ^ + H(t) = ±oo, H is non-integrable 
on [0, t s ], t s arbitrary and thus this model universe implements exactly Condition SI of 
Theorem 13.2.11 

Another basic feature of this solution is that the second derivative of the scale factor 

is 

d(t) = q(q - l)Bti- 2 - — L — -j— -» -oo, as t - t s , (AAA) 

so that \E\, and consequently p, diverges at t s as it is expected for a sudden singularity. 
Note that p, s > for q 2 t~ 2 (a s — l) 2 > —k, and since (// + 3p) = —6a/a (see Eq. ( 12.1.51) ). we 
can have ji > and /x + 3p > throughout the evolution of such an expanding universe. 
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We can further prove that under the additional assumption p > 0, a universe exhibiting 
a sudden singularity is geodesically incomplete [67]. Firstly, we know (cf. [I], pp. 141-2) 
that in an FRW universe with p > and p > 0, given any vector X at any point q, the 
geodesic 7(f) which passes through the point q = 7(0) in the direction of X is such that 

• either j(v) can be extended to arbitrary values of its affine parameter v, or 

• there is some value vq > such that the scalar invariant (R a p — l/2Rg a p)(R al3 — 
l/2Rg a P) is unbounded on 7([0,i>]). 

For a RW spacetime, the invariant (R a p — l/2Rg a p)(R a/3 — l/2Rg al3 ) is found to be 

12k 24a 2 k 12a 4 12a 12a 2 a 12a 2 

— + ^^ + ^ + — + ^^ + ^, 4.4.5 

a 4 a 4 a 4 a 6 a A a 2 

and since a — >■ a(t s ), H(t) — > — >■ 00, a/a — > —00 as £ — > £ s , we see that (i? a/ g — 
l/2Rg a/3 )(R al3 — l/2Rg a P) is unbounded at t s . Hence, we find that such a spacetime 
is geodesically incomplete. Note that this holds true irrespectively of the exact form of 
the scale factor, we only use the qualitative properties that define a sudden singularity 
(Definition EH]). 

Further, the relative asymptotic strength of a sudden singularity is immediately implied 
by the definition of this type of singularity and it is necessarily of the form (H ~ \D\ ~ 
a) « \E\. 

4.4.2 Fluids with p and /i connected 

Although Barrow in [5B] kept the pressure and the density unconnected in his example 
of a sudden singularity, Nojiri et al. in [52] focused on finding specific equations of state 
appropriate for the emergence of such a structure in a flat universe. An example of such 
an equation of state is given in [52] , 



p + p = -C(p - A*) 



(4.4.6) 
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where 7 is a positive constant. We are interested below in the asymptotic behaviour of 
the resulting solution of the cosmological equations in the limit of fi — > fi . Integrating the 
continuity equation, we find that the scale factor, given by 



(ao is an integration constant), remains finite as fi — > fio- From the Friedmann equation 
3H 2 = fi, it follows that H and therefore \D\ must also remain finite. However, in the 
limit fi — > /io, the pressure, and hence \E\ blow up. Solving the Friedmann equation for 
H, substituting in the continuity equation fi = 3HC(fio — /i)~ 7 , and integrating we find 



which remains finite (since the integrand decreases to zero) even when fi — > /io (recall 
that 7 > 0). This fact leads to the conclusion that singular behaviour is met within a 
finite interval of proper time t and, as the following theorem shows, this type of behaviour 
continues to be valid also in curved models. 

Theorem 4.4.1 A necessary and sufficient condition for a sudden singularity of type 
(S3, N 2 , -B3) to occur at t in a universe filled with a fluid with equation of state p + fi = 
—C(fi — fi)~"' , 7 > ; is that fi — > fi at t . 

Proof. Integrating the continuity equation fi — 3CH(fi — /i)~ 7 = 0, we find 




(4.4.7) 






(4.4.9) 



which is finite as t — > t . Also since 




(4.4.10) 
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\E\ 2 = 3- 2 =^ + 3p) 2 = ±-(2^ + 3C(^-^) 2 , (4.4.11) 
cr 12 12 

\D\ 2 = 3 (V + 4Y = C (4-4.12) 



we see that as t — > to, H and \D\ remain finite whereas \E\ diverges, leading to a sudden 
singularity. The converse is immediate. ■ 



4.4.3 Milder-than-sudden singularities 

A sudden singularity is a very mild type of singularity in the sense that the only diverging 
physical quantity that characterises it is the \E\ component of the Bel- Robinson energy. 
The only other milder type of singularity that could occur would be one with a bounded 
Bel-Robinson energy as this is described by the B± case in our classification scheme. This 
can indeed happen. Consider a flat universe filled with a fluid with equation of state (cf. 
[52]) 

ABu 213 ^ 1 

P + l i = - AfJ p-i + B > °</ 3<1 /2. (4.4.13) 

Such a fluid, admits the exact solution a = a e T8/3 when (3 = 1/5. Then H = (8/3)r 5 / 3 , 
H = (40/9)r 2 / 3 and H = (80/27)r~ 1 / 3 , anc i t h us as T _^ q, a, a, a, H, H all remain finite 
whereas H becomes divergent. We can see that in this model the Bel Robinson energy at 
the initial time, 13(0), is finite whereas its time derivative is given by 



B(t) 



d , ■■ „ ■ . / k „ \ / kH 



(4.4.14) 



hence B(r) — > oo at r — > 0. Since the derivative of the Bel-Robinson energy diverges we 
may interpret this singularity geometrically as one in the velocities of the Bianchi (frame) 
field. At t s , the Bianchi field encounters a cusp and its velocity diverges there |70J. We 
believe this to be the mildest type of singularity known in the literature to date. 

Another example of a very mild singularity is given by Borde et al in [89]. These 
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authors define an averaged-out Hubble function, H av (given by Eq. (14.4.181) below), along 
a null or timelike geodesic and showed that every past-directed null or timelike geodesic 
satisfying H av > throughout the past is not past complete. To see this consider first a null 
geodesic with affine parameter A. Then using Eq. (I2.2.5P and substituting the Christoffel 
symbols (see footnote ([1]) in p. 14), we find that the geodesic equation is given by 

d 2 t a f 9 / dr\ 2 „ 9 / d9\ 2 



+ ~ \ a TT + f a TT +a 2 fsm 2 6 ^ =0. (4.4.15) 



d\ 2 a \ \dX J \dX J \dX / 

Using the fact that the geodesic is null (ds 2 = 0) we have that 



d 2 t a f dt\ 2 , 

+ -[—) =0 (4.4.16) 



and by integration we find that 



d\ 2 a \dX J 



dX oc a{t)dt. (4.4.17) 



Setting dX = (a(t)/a(t s ))dt, the affine parameter can be normalised so that dX/dt = 1 for 
t = t s .The averaged-out Hubble function is defined by 

Now by considering those null geodesies that satisfy H av > 0, we find 

1 f x{u) TT/XN „ 1 f a(ts) da 1 

< H av = T771 V7TT / H(X)dX = — — — / — - < 



X(t s ) - \(U) J KU) X(t s ) - \{U) J a{u) a(t s ) X(t s ) - \(U) ' 

(4.4.19) 

This shows that the affine parameter of every past-directed null geodesic having H av > 
must take values only in a finite interval, which implies geodesic incompleteness. A 
similar proof is obtained for the case of a timelike geodesic. Hence, all causal past-directed 
geodesies obeying the condition H av > are past incomplete. 

Note that for null geodesies, Eq. (14.4. 19[) implies that H is finite but only in a finite 
time interval (an analogous conclusion holds for a timelike geodesic) and this is no other 
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than the S3 (necessary) condition for a singularity of Theorem 13.2. II We therefore conclude 
that under the assumptions of [89J, privileged observers cannot exist in such expanding 
universes for an infinite proper time because had they existed Theorem 12.2.21 would imply 
that spacetime ought to be timelike and null geodesically complete. 



4.5 Sudden-rip singularities 

A sudden-rip singularity arises once we allow not only \E\ but also \D\ to diverge at a 
finite time, while the scale factor remains finite there. This of course implies that H, p 
and fi become also divergent. This possibility was explored in some examples in [52]. An 
appropriate equation of state to illustrate this type of asymptotic behaviour is to take [52] 

p + H=-BfJ , 1 (3>1. (4.5.1) 

Integration of the continuity equation fi — ?>BHji^ = then gives directly [52] 

a = a exp(^-^), (4.5.2) 

which is finite everywhere since we have chosen /3 > 1. We can find the relation between 
the proper time and the density by first substituting Eq. (14.5.11) in the continuity equation 
and then using the Friedmann equation 3if 2 = /i. We obtain 



or equivalently [52J 

2 u- p+l12 

t ~t + 7T- (4.5.4) 

Therefore, when fi — > 00, t — *■ to and so \p\, H, \D\ and \E\ all diverge within the finite 
time to, meaning that this is an (S2, N 2 , B\) singularity. The asymptotic strength of this 
singularity is found to be: a « H « \D\ « \E\. 
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This behavior is a special case of the following general result which holds also in curved 
models. 

Theorem 4.5.1 A necessary and sufficient condition for a sudden-rip singularity (of type 
(5*2, N 2 , Bi) ) to occur at t s in an isotropic universe filled with a fluid with equation of state 
p + ji = —Bfi 13 , (3 > 1, is that \i — > oo as t — > t s . 

Proof. From the continuity equation we have 

a = a ° 6XP ( 3 £(1 -/?) )' (45 - 5) 

and so a — > a as t — > £ s . Since 

# 2 = ^-^ l^| 2 = ^(2/x + 35/) 2 , |£| 2 = ^, (4.5.6) 
6 a z ±2 3 

we see that as t — > t s , H, \E\ and \D\ diverge leading precisely to a sudden- rip singularity. 
The converse is obvious. ■ 



4.6 Discussion 

We have illustrated in this chapter that our methods of classifying singularities through the 
behaviour of the Hubble parameter the scale factor and the Bel- Robinson energy provide a 
rigorous basis for a useful description of the dynamics of a variety of cosmological models 
that exhibit finite-time singularities. Even though these model universes come from com- 
pletely different motivations yet they reveal their tendencies to have finite-time singularities 
in the sense of geodesic incompleteness for primarily one reason, the non-integrability of 
H. Finite-time singularities such as collapse, big-rip or sudden singularities occur because 
H fails to be integrable for an infinite time interval. In such universes, privileged observers 
cannot exist for an infinite proper time. We know from Theorem 12.2 . 21 that if such observers 
existed for an infinite proper time, then the universe would be timelike and null geodesi- 
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cally complete. Indeed, this is the case in various families of complete models ranging from 
dark energy, to inflationary to braneworlds that we have already studied. 

Naturally, the fact that the non-integrability condition for H can be sustained in differ- 
ent models stems from the various ways that these models are constructed and the specific 
features they share. This non-integrability may be supported either by the presence of 
phantom-dark energy driving the observed cosmic acceleration and leading to a divergence 
of both \E\ and \D\, or by a fluid with only a diverging \E\ as in [56], or in the inflationary 
character of the specific model as in [89], etc. All these different particularities are acknowl- 
edged and brought to light with the use of the Bel-Robinson energy. We have shown that 
with the combination of our classifying methods we have succeeded to obtain a complete 
picture of the reasons underlying a wide variety of possible dynamical behaviours of the 
universe on approach to a spacetime singularity. 



Chapter 5 



Braneworld singularities 

5.1 Introduction 

In previous chapters we studied singularities that arise within the framework of a four- 
dimensional RW-spacetime that models the whole entity of the observed universe. In this 
chapter we study possible singularities that develop when we consider the different case 
in which the observed universe is modeled by a four-dimensional hypersurface situated at 
a fixed position of an extra spatial dimension. The whole spacetime is five dimensional: 
there are four dimensions of space out of which only three are spanned by the hypersurface 
and one dimension of time; such a hypersurface is called a 3-brane while the full higher 
dimensional space is called the bulk. 

Although a brane can trap particles and forces making it impossible for them to escape, 
it does interact with the matter and forces in the bulk. This happens because gravity 
naturally extends over all dimensions of space and time and it can therefore influence the 
fields that are constrained to live within the boundaries of the brane [139] . Other matter 
fields in the bulk can also interact with the fields on the brane, and the strength of this 
interaction is controlled by a coupling function that is mo del- dependent [139J. 

Branes were discovered within the framework of string theory as locations where open 
strings are constrained to end. Their importance was reinforced when it was realised that 
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they can play a catalytic role in proving duality between the different versions of string 
theory. Driven by these new discoveries, a great deal of recent research focuses on building 
model universes that include branes (see Refs. [140j - [T9"T] for a partial list). 

All these different model universes that flourish in the recent literature are termed 
braneworlds and in many cases provided satisfactory explanations to some of the most 
challenging unresolved issues of cosmology and particle physics such as the hierarchy prob- 
lem and the cosmological constant problem. Below, we refer briefly to the ideas behind 
three of the most recognised braneworld cosmological scenarios and to the implications 
they offer for the singularity problem. The first scenario is the one introduced by Nima 
Arkani-Hamed, Savas Dimopoulos and Gia Dvali in |198j dubbed as the ADD scenario and 
the two scenarios that follow are those introduced by Lisa Randall and Raman Sundrum 
in [199] and [200] known as RS 1 and RS 2 . 

According to the ADD scenario all standard particles and forces except gravity are 
confined on a brane while the extra dimension is as big as a tenth of a millimeter; this is 
surprising since until then it was believed that extra dimensions were significantly smaller. 
In this scenario the hierarchy problem, that poses the question of why gravity at small 
distances appears to be so much weaker than the other forces, is solved: the effect of 
gravity on particles could fade away because of the large extra dimension so that gravity 
would appear much weaker than the brane-restricted electromagnetic strong and weak 
forces whose strength remains unaffected exactly because they cannot experience the extra 
dimension. 

A different approach to the hierarchy problem is provided in [199j and [200] . In the 
RS 1 scenario the extra dimension is considered to be small in size and bounded between 
two branes called the Gravitybrane and the Weakbrane. In this model the authors make 
essential use of a warped geometry where the warp factor that defines how space, time, 
energies and masses change with distance (measured along the fifth extra dimension) de- 
creases exponentially along the fifth dimension. A consequence of this configuration is that 
the energy of the bulk and the branes warps the spacetime in such a way that gravity is felt 
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strong on the Gravitybrane and weak on the Weakbrane which accommodates all standard 
particles and forces. In contrast, in the RS 2 scenario where there is only one brane, the 
extra dimension is infinite in size but it can nonetheless remain invisible provided that the 
spacetime is suitably warped. 

In this chapter we use the concept of a warped geometry in our braneworld models 
but we will not focus so much on the exploration of the phenomenological implications 
that these models offer, rather on the interesting features of finite-distance singularities 
that they exhibit as well as on very basic questions (to which we refer below) that the 
development of such singularities erases. 

In particular, we will focus on the analysis of a new version of a braneworld model 
that was originally proposed by Arkani-Hamed et al. in [73]. In the original version [71] 
the model consists of a flat three-brane embedded in a five dimensional bulk space that 
is filled with a scalar field minimally coupled to the bulk and conformally coupled to the 
fields on the brane. The metric describing the five dimensional spacetime is a warped 
product as in the RS scenarios. A basic feature of this model is the development of a 
singularity at a finite distance from the brane around which the warp factor, its derivative 
and the energy density of the scalar field all diverge. This singularity was used in [71] in 
order to provide an explanation to the cosmological constant problem while preserving four 
dimensional flatness. In particular, the authors in [71] argue that this singularity on the 
one hand successfully cuts off the size of the extra dimension giving rise to the observed 
four dimensional flatness, and on the other hand it acts as a reservoir through which all of 
the vacuum energy is emptied so that its value is brought to agrement with the observed 
one. 

Our interest in this particular braneworld model springs precisely from the fact that it 
exhibits a finite-distance singularity which bares resemblance to the finite-time singularities 
that are the core of our research. Our quest is to provide answers to the following questions 
that naturally arise: How general is the behaviour found in [H]? Is the finite-distance 
singularity the only possible outcome when we pass on to a curved instead of a flat brane? 
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How strongly are the dynamical features of the model affected by a different form of bulk 
matter such as a perfect fluid? Is the type of singularity found in [71] the only possible one 
in these new braneworld models, and if not, what other types of finite-distance singularities 
can arise? 

In this chapter we show that in a curved brane the finite-distance singularity found 
in [71] shifts to an infinite distance away from the brane, while with a perfect fluid in the 
bulk new types of finite-distance singularities emerge. On approach to these singularities 
the warp factor can do only two things: It will either tend to zero (collapse-brane I or 
II), or diverge (big rip-brane singularities^) while the energy density of the bulk-matter 
components usually diverges. Our analysis implies that the effect of an infinite energy leak 
to the extra dimension is indeed generic. For a qualitative analysis of such models the 
method of asymptotic splittings, expounded in Appendix E, is a most useful tool. 

The plan of this chapter is as follows: In Section 5.2, we review the basic aspects 
of the embedding of a three-brane in a five dimensional bulk space. In Section 5.3, we 
give the field equations that describe the dynamics of a three-brane embedded in a five 
dimensional bulk with a scalar field; in Subsection 5.3.1, we analyse the possible asymptotic 
behaviours of that model using the method of asymptotic splittings. In Section 5.4, we 
explore the consequences of altering the matter in the bulk. In particular, the case where 
the bulk is filled with a perfect fluid has a dramatic effect on the global dynamics. We no 
longer encounter just the singularities of the type found in [71] but now we have dynamics 
characterised by the finite-distance singularities referred to in the previous paragraph. In 
Section 5.6, we conclude and discuss our results. 



5.2 Brane geometry 

Consider a four- dimensional hypersurface S that is situated at the arbitrary but fixed 
position Y* where Y is the coordinate along the fifth dimension of the embedding space. 
The brane has negligible thickness compared to the extra dimension (and also the other 
^ee Def. I5XT1 
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dimensions) extending over the higher dimensional bulk space and can be viewed as a 
thin surface (surface layer). The brane therefore naturally divides the bulk in two regions 
V + and V~ each containing one of the boundary surfaces S + and S~ that confine the 
brane. The two regions V + and V~ are expected to look identical from S + and X - , 
the only difference being that their normal vector should be reversed; in other words our 
braneworld admits aF-* — Y symmetry. The extrinsic curvatures of the two embeddings, 
K^ b and K~ b , therefore satisfy 

K^ b = —K~ b , everywhere on X. (5.2.1) 

Using the Gauss-Codazzi formalism, analysed in Appendix B.2, we can write the Ricci 
5-tensor in a form analogous to ( IB. 2. 201) (note that e{n) = 1 ), i.e., 

R ab = (4) Rab ~ K ab K + 2K an Kl l - (5.2.2) 

or 

Rab = + Z ab , (5.2.3) 

where 

Zab = (4) R ab ~ K ab K + 2K an K?. (5.2.4) 
Tracing the Einstein equations 

Rab - ^g a bR = K 2 5 T ab , (5.2.5) 

where k\ = 871G5 = M^ 3 with G5 being the five dimensional Newton's constant and M5 
the five dimensional Planck mass, we find that the scalar curvature is given by 

R=- 2 -4T. (5.2.6) 



Substituting back in (I5.2.5P leads to the following equivalent form of the Einstein equations: 
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Rab — «5 ( Tab — -^QabT 



(5.2.7) 



Inserting Eq. (I5.2.3P in Eq. (I5.2.7P and integrating over the layer defined by the thickness 
of the brane we obtain [201], [202] 



JY*-e V ° / 

In the passage to the limit e — > and for fixed —K^ h) K~ b we naturally expect that 
K a b remains bounded inside the layer and hence the integral of Z a b over the layer should 
vanish [201] . [202] . Taking thus the limit of the above equation for e — > and using Eq. 
(15.2.11) . we find the following junction condition [201] . [202] (see also [203] pp. 551-556) 



that relates the jump of the extrinsic curvature [K a b] across the brane to the surface-energy 
momentum tensor Sab (this follows from the discussion below) defined by 

/■Y*+e 

S ab = lim / T ab dY, (5.2.10) 

with S = g ab S a b being the trace of S a b- It can be proved (cf. [201] . [202] . and [203] pp. 
551-556) that integration of the rest of the components of the Einstein equations gives 
S55 and S$ a with a = 1,2,3,4 both zero and and hence, we conclude that there is no 
flow of energy from the fields restricted on E into the bulk. It is thus implied that the 
quantity S a b acts as a delta function source of matter existing only on £ and vanishing 
off S. This last property suggests that S a b plays the role of a surface-energy momentum 
tensor. Combining this with Eq. (15.2.91) . it follows that the extrinsic curvature of the brane 
is completely determined by the matter confined inside the brane [204] . 




(5.2.8) 




(5.2.9) 
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5.3 Asymptotics of scalar field-brane configuration 



Consider the braneworld model consisting of the three-brane configuration analysed in the 
previous Section and filled with a scalar field minimally coupled to the five-dimensional 
bulk and conformally coupled to any fields restricted on the brane. The total action Stated 
then splits in two parts, namely, the bulk action Sbulk and the brane action Sbrane, i-e., 



S total — Sbulk + Sbrane-, (5.3.1) 



with 



Sbulk = j d A xdY^~ 5 f A - ^(V0) 2 ) , (5.3.2) 
Sbrane = - J <fixy/gif{(j>), at Y = Y m , (5.3.3) 

where A is a parameter defining the type of scalar field (ft, n\ = M^ 3 , M 5 is the five 
dimensional Planck mass, and f(4>) denotes the tension of the brane depending on the 
scalar field. The tension of the brane can be understood physically as the resistance that 
one encounters when trying to distort the brane: zero tension means that the brane has 
no resistance so that any distortion would influence its form dramatically, while infinite 
tension would deprive the brane's liberty to evolve dynamically [139] . 

Varying the total action (15.3.11) with respect to g AB , we find the five- dimensional 
Einstein field equations in the form 



Rab - \ 9 abR = X4 (va0V b - \qabWT) - ^^&* S <X)$ ( 5 - 3 - 4 ) 
while the scalar field equation is obtained by variation of the action ([5.3. ip with respect to 



2 Note here that by working in a similar way as in Appendix D, we may interpret the energy-momentum 
tensor of the scalar field as that of a perfect fluid in the bulk having density p = p(Y) and pressure 
P = P(Y), i.e., T A b = {p + P)u A u B - Pqab, where u A = (0,0,0,0,1), so that P = 1/2 (V0) 2 = p. 

3 The left-hand side of this equation is the analogous left-hand side of the equation of motion (1D.0.7|) 
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1 SSf) rane 



8{Y), 



(5.3.5) 



where A, B = 1,2,3,4,5 and a,f3 = 1,2,3,4 while S(Y) — 1 at Y — Y* and vanishing 
everywhere else. 

We assume a bulk metric of the form 



ds 2 = a 2 (Y)ds 2 + dY 



(5.3.6) 



where ds 2 is the four dimensional flat, de Sitter or anti de Sitter metric, i.e. 



ds 



7.2 



-dt 2 + dr 2 + r 2 da 2 



-dt 2 + 



cosh 2 (Hi) 

m 

cos 2 (Ht) 



(dr 2 + sin 2 rda 2 



dt 2 + co^Ht)^ dr 2 +sinh 2 rda 2^ 



(5.3.7) 



where da 2 = d6 2 + sm 2 9d(p 2 . The field equations then take the form 



a 
a 



2 A'2 



kH 2 



12 a 2 



0. 



(5.3.8) 
(5.3.9) 
(5.3.10) 



where k = 0, ±1 and H' 1 is the de Sitter curvature radius and the prime denotes dif- 
ferentiation with respect to Y. The variables to be determined are a, a' and (p. These 
three equations are not independent since Eq. (15.3.91) was derived after substitution of Eq. 
(15.3.81) in the field equation G aa = K^T aa , a = 1, 2, 3, 4, 



kH 2 



a a 
a a 2 a? 



','2 



-Xnir—-. 



(5.3.11) 



In our analysis we use the independent equations (I5.3.9P and (15.3.101) to determine the 
of a scalar field in a four-dimensional spacetime with V(4>) — (see appendix D). 
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unknown variables a, a' and <$' and the Eq. (15.3.81) will then play the role of a constraint 
equation for our system. 

The extrinsic curvature of the three-brane is calculated to be 

K a p = \^9ap = a(Y)a(Y), (5.3.12) 

while the surface energy-momentum tensor (defined only on the brane and vanishing off 
the brane) is now given by 

s <# = -*%H£r WW = -*<X.)fW.)). (5-3.13) 

v 95 "9 



Solving Eq. (15.3.41) on the brane is equivalent to substituting Eqs. (15.3.121) and (15.3. 13j) in 
condition (15.2.91) which leads to the following relation to be satisfied by a'(Y*): 

a'%) = _^/(0(y.)) o (n). (5-3.14) 
o 

On the other hand, Eq. (I5.3.5P gives directly 

A ° 'VtoV) = -U f dx 4 Jgtf'(<j>)6(y). (5.3.15) 
^/9lJ 



Integrating over the interval [K* — e, Y* + e] and taking the limit as e — > 0, we finally find 

W) = (5.3.16) 
where we have used the Y — > — Y symmetry, i.e., /+ = (j)'(Y* + e) = —(f)' (Y* — e) = — <j)'~ . 

Our purpose now is to find all possible asymptotic behaviours around the assumed 
position of a singularity, denoted by Y s , emerging from general or particular solutions of the 
system of equations ( 15.3.8l) -( !5.3.10p describing our braneworld cosmology. As we mentioned 
in the Introduction a most useful tool for this analysis is the method of asymptotic splittings 
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([76], see also Appendix E). According to this method, we start by setting 

x = a, y = a , z = <fi', (5.3.17) 

where the differentiation is considered with respect to T = Y — Y s . The field equations 
(15.3.91) and (15.3.101) become the following system of ordinary differential equations: 

x' = y (5.3.18) 
y' = -\Az 2 x (5.3.19) 

z' = -Ay-, (5.3.20) 

x 

where A = k\/A. Hence, we have the vector field 

f = (y, -\Az 2 x, -%^) T • (5.3.21) 

Equation (15.3.81) does not include any terms containing derivatives with respect to T; it is 
to be regarded as a constraint equation for the above system (I5.3.18p - (l5.3.20p . In terms of 
the new variables the constraint is written in the form 



£ = + (5-3.22) 



Substituting the forms 



where 



(x,y,z) = (aT p ,(3T^ST r 



(p, q, r) E Q 3 and (a, [3, 5) E C 3 \ {0}, (5.3.23) 

in the dynamical system (I5.3.18p - (l5.3.20p . we seek to determine the possible dominant 
balances in the neighborhood of the singularity, that is pairs of the form 



£> = {a, p}, where a. = (a, ft, 5), p=(p,q,r). 



(5.3.24) 
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For our system fl5.:-i.l8D - fl5.:i2()l) we find*: 



B 1 = {(a, a/4, V3/(4>/AA)), (1/4, -3/4,-1)} 
B 2 = {(a, a, 0), (1,0,-1)} 



(5.3.25) 
(5.3.26) 



Since (15.3. 18[) - (I5.3.20j) is a weight-homogeneous system, the scale invariant solutions 
given by the above balances are exact solutions of the system. The balance B\ satisfies the 
constraint Eq. (I5.3.22p only for k = corresponding thus to a general (as we will see in the 
next paragraphs) solution that allows only for a flat spatial brane geometry. In contrast, 
B 2 corresponds to a particular solution for a curved brane since it satisfies Eq. (15.3.221) for 
k 7^ and a 2 = kH 2 (here we have to sacrifice one arbitrary constant by setting it equal 
to kH 2 ). 

Let us now focus on each of the two possible balances separately and build series ex- 
pansions in the neighborhood of the singularity. We start by calculating the Kowalevskaya 
exponents, eigenvalues of the matrix /C = -Df(a) — diag(p), where -Df(a) is the Jacobian 
matrix of /, 

( 1 \ 



Df(x,y,z) 



V 



-AXz 2 
Ayz 







-2A\xz 
_4y 

x 



(5.3.27) 



/ 



ar x 

evaluated on a. For the first balance we have that a = (a, a/4, y/3/ (Ay/ AX)), and p 
(1/4,-3/4,-1), thus 



/ 



K 



1 

~4 
3 

~W 
V3 



1 

3 
4 



\ Aa\/AX ay/AX 







V3AX 



a 



(5.3.28) 



3 ^ We have ignored here a third balance B3 = {(a, 0, 0), (0, — 1, —1)} since it does not give the necessary 
— 1 as one of the Kowalevskaya exponents (see Appendix E for more details). 
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The /C-exponents are then given by 

spec(/C) = {-1, 0, 3/2}. (5.3.29) 

As it follows from the analysis in Appendix E, these exponents correspond to the indices of 
the series coefficients where arbitrary constants first appear. The —1 exponent signals the 
arbitrary position of the singularity, Y s . This first balance has two non- negative rational 
eigenvalues which means that it describes a general solution, i.e., a solution having the 
full number of arbitrary constants. In order to construct an asymptotic expansion of this 
solution, we substitute in the system fl5.3.18l) - fl5.3.20l) the series expansions 

x = T p (a + £°l lCj T J >), 

where x = (x,y,z), Cj = (cji, Cj2, C73), and s is the least common multiple of the denomi- 
nators of the positive eigenvalues; here s — 2, and the corresponding series expansions are 
given by, 



x 



£,°lo9i T ' /2+1/4 > V = £,°Vi2T j/2 - 3/4 , z = ^f =0 c j3 T^ 2 -\ (5.3.30) 



Following this, we arrive at the asymptotic solution around the singularity in the form 

x = aT 1 / 4 + ^c 32 T 7 / 4 + • ■ ■ (5.3.31) 
^T~ 3 / 4 + c 32 T 3 / 4 + ■ • • (5.3.32) 

^ :T -i _ ^L C32T i/2 + . . . . (5.3.33) 



y 



AVAX 7 a VAX 

The last step is to check if for each j satisfying j/2 = p with p a positive eigenvalue, 
the corresponding eigenvector v of the K matrix is such that the compatibility conditions 
hold, namely, 

v T ■ Pj = 0, (5.3.34) 
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where Pj are polynomials in Cj, . . . , c,_i given by 



K,Cj - (]/s)cj = Pj. 



(5.3.35) 



Here the corresponding relation j/2 = 3/2, is valid only for j = 3 and the associated 
eigenvector is 



aVAX 7aVA\ 



v 



Vs ' V3 ' 

The compatibility condition, v • {K, — (3/2)X 3 )c 3 = 0, therefore indeed holds since 



(5.3.36) 



/ 



{K - (3/2)X 3 )c 3 = c 32 



7 
~4 

3 
v/3 



3 

>/3 



a 







-\73AA 



V Aa\fA\ aVAX 



3 
2 



4 
7 



4y/3 

V 7av / AA / 



V / 

(5.3.37) 



This shows that a representation of the solution asymptotically with a Puiseux series as 
this is given by Eqs. fl5.3.31l) - fl5.3.33l) is valid. We thus conclude that near the singularity 
at finite distance Y s from the brane the asymptotic forms of the variables are 



a — > 0, a — > oo, 



oo. 



(5.3.38) 



This is exactly the asymptotic behaviour of the solution found previously by Arkani- 
Hammed et al in [74]. 

However, the previous behaviour is not the only possible one; the analysis below shows 
that a very different outcome results by the asymptotic analysis of the second balance. 
The /C-matrix of the second balance is 
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K = Df((a, a, 0)) - diag(l, 0, -1) 



-1 1 

0-3 



(5.3.39) 



with eigenvalues 



spec(/C) = {-1,0,-3}. 



(5.3.40) 



For the balance £>2 we find two distinct negative integer /C-exponents and and hence we 
expect to find an infinite negative expansion of a particular solution (recall that we had to 
sacrifice one arbitrary constant) around the presumed singularity at Y s with the negative 
/C-exponents signaling the positions where the arbitrary constants first appear [205J. We 
therefore expand the variables in series with descending powers of T in order to meet the 
two arbitrary constants occurring for j = —1 and j = —3, i.e., 



(5.3.41) 



Substituting these series expansions back in the system ( 15.3.18l) -( 15.3.20l) . we find the fol- 
lowing asymptotic behaviour 



x = aT + c_xi + 
y = a H 

z 



C-33T- 4 + 



(5.3.42) 
(5.3.43) 
(5.3.44) 



We now check the compatibility condition for j = —3 (the associated eigenvector here is 
v = (0,0,1)). We find that 



P_ 3 = (/C + 3J 3 )c_c 



I 2 1 \ / ^ 

3 




/ 



V C33 / 






V / 



(5.3.45) 
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so that the compatibility condition is indeed satisfied. The expansions given by Eqs. 
fl5.3.42p - fl5.3.44p are therefore valid and we can say that as T — > 0, or equivalently as 
S = 1/T — > oo, we have that 

a — > oo, a' — > oo, </>'—> oo. (5.3.46) 

We thus conclude that there exist two possible outcomes given by Eqs. (I5.3.38[) and (I5.3.46j) 
for these braneworld models, the dynamical behaviours of which strongly depend on the 
spatial geometry of the brane. For a flat brane the model experiences a finite-distance 
singularity through which all the vacuum energy decays (since 0' — > oo, as 7 — > Y s ), 
whereas for a de Sitter or anti-de Sitter brane the singularity is now located at an infinite 
distance. Flat solutions correspond to the particular coupling used by Arkani-Hammed et 
al in [74J, the only choice to make this possible. This easily follows by using equations 
(15.3. 14p and (15.3.161) and solving the FRW equation (15.3. 8p on the brane for kH 2 : 

Clearly, k is identically zero if and only if: 




or equivalently, if and only if /(</>) oc e 2 ^^ K5 ^ (Arkani-Hammed et al in [71] have A = 3, 
and hence the appropriate choice for the brane tension is f(4>) oc e 2K5 ^). 

By working with other couplings we can allow for non-flat, maximally symmetric solu- 
tions to exist and avoid having the singularity at a finite distance away from the position 
of the brane. 

In the following Section we replace the scalar field in the bulk with a perfect fluid and 
explore the outcome of the resulting dynamical system. 
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5.4 Perfect fluid in the bulk 

We consider here the case that the bulk matter takes the form of a perfect fluid with 
equation of state P = 7p, where P and p are functions of Y only. This fluid can be viewed 
as an entity in its own right but we can also use it to represent an average contribution of 
the various scalar fields living in the bulk. 

The field equations describing the braneworld model in this case are derived by replac- 
ing the energy-momentum tensor in the first term of the RHS of Eq. (I5.3.4p by that of a 
perfect fluid, Tab = (p + P)uaUb — PgAB, ua = (0, 0, 0, 0, 1). They read: 

a - - (5.4.1) 
a o 

p' + 4(1 + i)Hp = 0, (5.4.2) 

— = -Mp + — . 5.4.3 

cr o cr 

Introducing the new variables 

x — a, y = a , w = p, (5.4.4) 

Eqs. f l5.4.ip - fl5.4.3|) can be written in the form 

x' = y, (5.4.5) 
y' = -2A ^ + ^ wx, (5.4.6) 
w ' = _4(i + 7 )^ W) (5.4.7) 



x 



with the constraint equation given by 



The system f l5.4.5l) - fl5.4.7p is weight homogeneous with vector field 
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y, -2A ^ 1 + 2 ^ wx, -4(1 + 7)-w 
3 x 

and by setting (x,y,w) = (aT p , (3T q , 5T S ), we find the following set of five possible bal- 
ances: 



7 



B 1 = { (a,ap,^-p 2 ),(p,p-l,-2)\, p = - 1 - , 7 ^ -1/2, -1,(5.4.9) 



2^;'^ ' y J' ' 2(7 + 1) 



7 £ 2 = {(a,a,0),(l,0,-2)}, 7 ^ -1/2, (5.4.10) 

-1/2B3 = {(a,a,0),(l,0,r)}, (5.4.11) 

-1/2B4 = {(«, a, 5), (1,0, -2)}, (5.4.12) 

_i /2 £ 5 = {(a,0,0),(0,-l,r)}, (5.4.13) 

where = 7= _i/ 2 £>j. These balances are exact solutions of the system and they should 

therefore also satisfy the constraint equation (15.4.81) . This fact alters the generality of the 
solution represented by each balance and determines the type of spatial geometry that it 
admits: The balances 1 B\ and _ 1/2*65 describe a potentially general solution corresponding 
to a flat brane, while the balances 7 i3 2 and -1/2-63 describe particular solutions of curved 
branes since we already have to sacrifice the arbitrary constant a by accepting that a 2 = 
kH 2 . The balance _i/ 2 £>4 on the other hand describes a particular solution of a curved or 
flat brane since we have to set 5 = (3/(2v4))(l — kH 2 /a 2 ) in order to satisfy Eq. (15.4.81) . 

Before we proceed to the analysis of each of these balances in the next Subsections, 
we find it useful to introduce the following definition of finite-distance singularities: 

Definition 5.4.1 Let a be the warp factor in the metric 115.3.6}) . a' its derivative and g 
the energy density of the matter component in the bulk. We call a finite- distance brane- 
singularity occurring with: 

i) a — > 0, a' —>■ 00 and g — > 00; a collapse-brane I singularity. 

ii) a — > 0, a' — > c, where c is an arbitrary constant, and g — > 00: a collapse-brane II 
singularity. 



82 



CHAPTER 5. BRANEWORLD SINGULARITIES 



Hi) a — > oo 7 a' — > oo and o — > oo: a fog rip-brane singularity. 

Note that with this terminology, the finite-distance singularity we met in the previous 
Section is a collapse-brane I singularity. 



5.4.1 The balance 7 B\ 

We begin our analysis by calculating the Kowalevskaya-matrix for the first balance. The 
Jacobian matrix of f is given by 



Df(x,y,z) 



( 



--(l + 2 1 )Aw 







\ 



--(l + 2 7 )Ar 



V 4(l + 7 )^ -4(1 + 7)- -4(1 + 7)- J 

X A X X ' 



(5.4.14) 



and 



y /Ci = Df[a, ap, —p 2 ) - diag(p, p - 1, -2) 



Df a, 



2(l + 7 )' 8A(l + 7 ) ; 



— diag 



1 l + 2 7 

2(1 + 7)' 2(1 + 7)' 



-2 



2(1 + 7) 



l + 2 7 



\ 4(1 + 1 ) 2 Aa 
The eigenvalues of this matrix are 



\ 



— \ + 21 \ --(1 + 2 7 )A* 

4(l+7) 2 2(1 + 7) 3 V [; 



2(1 + -f)Aa 



(5.4.15) 



/ 



spec( 7 /Ci) = < —1, 0, 



27 



1 + 7 



(5.4.16) 



5.4. PERFECT FLUID IN THE BULK 



83 



The last eigenvalue is a function of the 7 parameter and it is positive when either 7 < — 1 
or 7 > —1/2. We are going to focus on each of these two ranges of 7 and explore the 
possible outcome of the system and the character of the finite-distance singularities that 
occur each time. 

Consider first the case with 7 < — 1. For purposes of illustration, we take 7 = — 2. 
Then the balance and the _2/Ci-exponents read respectively, 



_ 2 £a = {(a,-a/2,3/(8A)),(-l/2,-3/2,-2)}, 
spec(_ 2 /Ci) = {-1,0,3}. 



(5.4.17) 
(5.4.18) 



Substituting the value 7 = — 2 in the system given by Eqs. fl5.4.5p - fl5.4.7l) and also the 
forms 



rp — yoo r . 1 T J ~ 1 / 2 11 — y°° r„P -3 / 2 7 
x — i^j =Q Cji 1 , y — i^j =0 Cj2 J- , z 



(5.4.19) 



we expect to meet the third arbitrary constant at j = 3. Indeed we find, 



x 



w 



8A 



T- 2 + c 33 T + 



(5.4.20) 
(5.4.21) 
(5.4.22) 



The compatibility condition is trivially satisfied for j = 3 since 



{K - 3J 3 )c 3 



2 

3 
4 



-- 2Aa 



c 3 3 



3 



5 , 

-Aa 





V / 



(5.4.23) 



V AAa 2Aa 

The series expansions given by Eqs. (15.4.201) - f l5.4.22f) are therefore valid asymptotically for 
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T — > so that 

a — ► oo, a' — > — oo, p — > oo. (5.4.24) 

We therefore conclude that the balance 7 £>i describes a general solution (corresponding to a 
flat brane) around a big rip-brane singularity when the perfect fluid satisfies a phantom-like 
equation of state, i.e., 7 < — 1. Note that using the analogy between the warp factor and 
the scale factor we can say that this singularity bares many similarities to the analogous 
one arising in the classical theory (see Section 4.3) since it is characterized by all quantities 
a, p and consequently P and H becoming divergent asymptotically. 

In order to illustrate the second range of values of 7, 7 > —1/2, we assume that 
7 = — 1/4. Then 

-1/4B1 = {(«,-2«/3,2/(3A)),(2/3,-l/3,-2)}, (5.4.25) 
spec(_ 1/4 /C 1 ) = {-1,0,2/3}. (5.4.26) 

Substituting in the system (I5.4.5|) - (|5.4.7p the particular value 7 = — 1/4 and the forms 

x = E- c il T^3+2/3 ) y = e« c,- 2 T^-V3 j z = E oc oCj . 3 t^ 3 - 2 , (5.4.27) 

we arrive at the following asymptotic expansions: 

x = a T 2 / 3 -^c 23 T 4 / 3 + --- , (5.4.28) 
V = ^T- 1 / 3 -^ac 23 T 1 / 3 + --- , (5.4.29) 
w = A T - 2 + C23 T- 4 / 3 + --- . (5.4.30) 



3A 



Again, the compatibility condition is bound to be satisfied since in this case 
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(JC - (2/3)J 3 )c 2 



4 
~3 


1 





2 


1 




~9 


~3 


3~ 


4 


2 


2 


3Aa 


Aa 


~3 



2 



C23 



2Aa 



\ 

Eqs. fl5.4.28l) - f!5.4.30p then imply that as T -> 0, 

a — > 0, a' — > cxd, p — >• 00. 



(5.4.31) 



; V / 



(5.4.32) 



This second asymptotic behaviour corresponds to a general solution of a flat brane that 
is valid around a collapse-brane I singularity. We thus regain a behaviour similar to that 
of the balance in (I5.3.25p . However, an important outcome of our analysis is that there is 
now in addition to that previously met, an equally general distinct big rip-brane singularity 
that may also arise during the evolution of a flat brane. 

5.4.2 The balance 7 £> 2 

For the second balance 7 £>2, we find that 



y /C 2 = Df (a, a, 0) - diag (1, 0, -2) 



and hence its spectrum is 



/-I 1 

\ 



\ 

|a«(1 + 2 7 ) 
-2(1 + 2 7 ) J 



(5.4.33) 



spec( 7 /C 2 ) = {-l,0,-2(l + 2 7 )}. 



(5.4.34) 



We note here that the third arbitrary constant will appear at j = 



-2(1 + 27). After 
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substituting the forms 



(5.4.35) 



in the system fl5.4.5l) - fl5.4.7p . we may try giving different values to 7: inserting 7 = —3/4 
in the system we meet the third arbitrary constant at j — 1 (spec(_3/4/C2) = { — 1,0, 1}). 
The variables (x, y, w) then expand in the following way: 



x = a 1 H C13 I 

6 

y = a + — Ci 3 T + - 
w = ci 3 T _1 H . 



(5.4.36) 

(5.4.37) 
(5.4.38) 



We ought to check the validity of the compatibility condition for j = 1. But this is again 
trivially satisfied since 



-3/4^2 — ^3)ci 



/ 



V 



-2 1 
-1 Aa/3 




N ^ ( Aa/6 ^ ^ n \ 



Cl3 



Aa/3 

V 1 J 






V / 



(5.4.39) 



The series expansions Eqs. f l5.4.36p - fl5.4.38p are therefore valid ones, and we conclude that 
as T — > 0, 

a — > 0, a' — > a, p — >• 00, (5.4.40) 

It will follow from the analysis below that the behaviour of p may vary according to our 
choice of 7. 

Indeed for 7 = — 1 (spec(_i/C 2 ) = {—1, 0, 2}) we find that 



x 

y 

w 



aT + — c 23 T 3 + 
9 

a + ~Y C 23± + ■ ■ 
C23 H • 



(5.4.41) 

(5.4.42) 
(5.4.43) 
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Note that the compatibility condition is satisfied also here: 

/ 



.i/C 2 - 2J 3 )c 2 



We see then that as T — > 0, 



-3 


1 





-2 









\ ( Aa/9 ^ 



C2 3 



Aa/3 

V 1 J 





V / 



(5.4.44) 



a — > 0, a' — ► a, p — > C23, a^O, 



(5.4.45) 



which shows that distinctly different behaviour than (15.4.401) . 

A yet different behaviour is met if we choose 7 = —5/4: The /C-exponents are in this 
case given by spec(_5/4/C2) = { — 1,0,3}, and the series expansions read: 



x 

y 

w 



Aa A 

aT + C33T 4 + ■ 

12 

a + —c 33 T 3 + --- 
C33T H . 



(5.4.46) 

(5.4.47) 
(5.4.48) 



These expansions are valid locally around the singularity since the compatibility condition 
holds true, 



/ 



-5/4^2 — 3X3) Ci 



For T — > we have that, 



-4 


1 








-3 


Aa 












\ 



C 3 2 



( Aa/ 12 ^ 

Aa/3 
1 



f \ 



V / 



(5.4.49) 



a — > 0, a — > a, p — > 0, a^0. 



(5.4.50) 



This balance describes therefore the asymptotic behaviour of a particular solution of 
a curved brane around a finite-distance singularity where a — > and a' — > a, i.e., a and a' 
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behave exactly as in the case of a collapse-brane II singularity but now the behaviour of the 
energy density of the perfect fluid varies dramatically from zero to infinity (p — ► 0,p s , oo) 
depending on the value of 7. This allows for new dynamical outcomes of the system 
f l5.4.5|) - fl5.4.8p : we can have a constant or even no flow of energy through the finite-distance 
singularity and into the extra dimension. 

A totally different result is obtained when we consider a family of fluids having a 7 
greater than —1/2. For 7 = 0, we have that spec(o/C 2 ) = {—1,0, —2} and hence we may 
expand (x, y, w) in descending powers in order to meet the arbitrary constant appearing 
at j — — 1 and j = —2, i.e., 



x 



Y~°°r .-,T'j +1 11 — Y-°°r n T3 
^5=0^1 1 1 V — 2 -'j=0 c j2 J- 



(5.4.51) 



We find: 



x 

y 

w 



Aa „_ 2 
a + — c_ 23 T H , 

c_ 23 T- 4 + -- - . 



(5.4.52) 

(5.4.53) 
(5.4.54) 



The compatibility condition is satisfied since 



)/C 2 + 2J 3 )c_ 2 



1 1 

2 -2Aa/3 




N ( -Aa/3 ^ / n \ 



C- 2 3 



v 1 y 






V / 



(5.4.55) 



AsS'=l/T— > 00 we conclude that 



a — > 00, a — > 00, p — » 00, 



(5.4.56) 



and we can therefore avoid the finite-distance singularity in this case. 
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5.4.3 The balances -1/2B3, -1/2&4 and -1/2B5 

In this Section, we analyse the behaviours that are described by the three balances which 
are valid once we set 7 = — 1/2 in the system fl5.4.5|) - fl5.4.7|) . We show that the balance 
_i/ 2 £> 3 , for r < —2, and the balance -1/2^4 imply the existence of a collapse-brane II 
singularity. However, setting r > —2 in the balance -1/2*83 leads to the avoidance of 
the finite-distance singularity. The balance -1/2B5 on the other hand, describes a smooth 
behaviour since the warp factor, its derivative and the energy density of the fluid, remain 
finite asymptotically. 

We begin with the analysis of the balance -1/2-63 . The K,- matrix for this balance is 



/ 



-1/2^3 = Di (a, a, 0) - diag(l, 0, r) 



-1 1 

-2-r 



\ 



(5.4.57) 



and hence, 



spec(_i/ 2 /C 3 ) = {-l,0,-2-r}. 



(5.4.58) 



Taking — 2 — r > 0, we have two non- negative /C-exponents. For r = — 3 we substitute the 
forms 



x 



y°° r-,T^ +1 11 — y°° r „p 7 — y°° r,T j_3 



(5.4.59) 



and arrive at the expansions 



x = aT + ■ ■ ■ , 

y = a H , 

w = ci 3 T~ 2 + - 



(5.4.60) 
(5.4.61) 
(5.4.62) 



The compatibility condition is satisfied since 
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\ 



-2 1 
0-10 




\ / n \ / n \ 

Cl3 






v 1 ; 






V / 



(5.4.63) 



The expansions ( 15.4. 67P -I 15. 4.69ft are therefore valid ones in the vicinity of the singularity 
and the general behaviour of the solution is characterised by 



a — > 0, a 1 — > a, p — > oo, 



(5.4.64) 



The balance -1/2B3 for r < —2 implies therefore the existence of a collapse-brane II singu- 
larity. 

If on the other hand we take r > —2, we have two negative /C-exponents. If we choose 
the value r = 0, then the spectrum is found to be 



spec(_ 1/2 /C 3 ) = {-1,0,-2}, 



(5.4.65) 



and so inserting the forms 



x 



y~°° r ., r ri+ 1 „, — y-00 yj _ y-00 r yj 



(5.4.66) 



we obtain 



x 

y 

w 



«T + c_i; 



a. 



(5.4.67) 
(5.4.68) 
(5.4.69) 



which validates the compatibility condition at j = 



—2 since 
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. 1/2 /C 3 + 2J 3 )c_ 2 



We see that as S = 1/T — > oo, 



/ 1 1 \ 

2 




C-2 3 





V 1 / 





V / 



(5.4.70) 



a— >c_n, a! —> a, p — >• oo, a^O, 



(5.4.71) 



so that the balance -1/2-63 for r > —2 offers the possibility of avoiding the finite-distance 
singularity met previously. 



The balance -1/2^4 is one with 



/ 



-V 



2/C4 = Df (a, a, 0) - diag(l, 0, -2) 



-1 1 





25 

a 







o\ 





■2 



and 



spec(„i /2 /C 4 ) = {-1,0,0}. 



(5.4.72) 



(5.4.73) 



We note that the double multiplicity of the zero eigenvalue reflects the fact that there were 
already two arbitrary constants, a and 5 in this balance (5 will have to be later sacrificed 
in order for that balance to satisfy the constraint Eq. (15.4.81) ). We can thus write 



x = aT + ■ ■ ■ 

V = a H , 

w = 5T- 2 + - 



(5.4.74) 
(5.4.75) 
(5.4.76) 



so that as T — > a collapse-brane II singularity arises, i.e., 



a — » 0, a' — > a, p — > 00, 



(5.4.77) 
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Finally we come to the balance -_i/ 2 £> 5 = {(a, 0, 0), (0, — 1, r)}. This has 



/ 



_i /2 /C 5 = £>f (a, 0, 0) - diag(0, -1, r) 



1 
1 

-r 



(5.4.78) 



and 



spec(_i /2 /C 5 ) = {1,0, -r}, 



(5.4.79) 



so we shall have to set r = 1 in order to have the —1 eigenvalue corresponding to the 
arbitrary position of the singularity, t s . After substitution of the forms 



x 



(5.4.80) 



we find that 



x = a + c\ i H , 

y = ciH , 

w = OH . 



(5.4.81) 
(5.4.82) 
(5.4.83) 



The compatibility condition is satisfied since 



(-1/2^5 - ^3) c l 



-1 1 





\ ( , \ 

Cll 



/ 



1 
1 

V / 



/o\ 



V / 



and we see that as T — > 0, 



(5.4.84) 



rr — > a, y — > Cn, w — > 0, «^0. 



(5.4.85) 



5.5. CONCLUSIONS 



93 



5.5 Conclusions 

The following two tables summarise the overall features of a braneworld consisting of a 
flat, or curved three-brane embedded in a five-dimensional bulk space with a scalar field 
(Table |5\TI) . or a perfect fluid (Table 1572]) . In what follows, we denote by f s the non-zero 
constant that is the limit of a function / on approach to the singularity at Y s , and by FDS, 
IDS a finite-distance and an infinite-distance singularity respectively. 

Table 5.1: Singularities in the case of a scalar field in the bulk 



Balance 


k 


Type of 


singularity 




Bi 


k = 


FDS: a - 


-> 0, a' — > oo, (ft 


— > oo 


B 2 


k^O 


IDS: a - 


-> oo, a' — > oo, 


— > oo 



Table 5.2: Singularities in the case of a perfect fluid in the bulk 



Balance 




7 


r 


Type of singularity 


•A 


fc = 


7 > -1/2 




FDS: a — > 0, a' — > oo, p — > oo 






7 < -1 




FDS: a — > oo, a' — > oo, p — > oo 


7 B 2 




7 < -1/2 




FDS: a — > 0, a' — > a, p — > 0, p s , oo 






7 >-l/2 




IDS: a — > oo, a' — ► oo, p — >■ oo 


-1/2^3 


k + 


7 = -1/2 


r < -2 
r > -2 


FDS: a — > 0, a' — > a, p — > oo 
IDS: a — ► oo, a' — > oo, p — >■ oo 


-1/2&4 


fc = 0,±1 


7 = -1/2 




FDS: a — >• 0, a' — > a, p — > oo 



For a flat brane and for both forms of bulk matter we find that there exists a collapse- 
brane I singularity as this is implied by the balances B\ and 7 £>i with 7 > —1/2. When 
we consider a scalar field this is the only possible type of FDS that may arise. In contrast, 
when a perfect fluid is present in the bulk there are in addition two new types of FDS: 
The first one is a very distinct big rip-brane singularity that appears only when a phantom 
equation of state is considered (7 < —1). This new possibility is offered by the balance 
y Bi and corresponds to a general solution. The second new type of FDS is that of a 
collapse-brane II singularity that is implied by the balance However this feature 

is less general than the ones described by the balance 7 £>i. This follows not only because 
we had to fix the value of 7 equal to —1/2 in the balance -1/2B4, but also because we had 
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to sacrifice one arbitrary constant by setting 8 = 3/ (2 A) in order to satisfy the constraint 
Eq. (15.4. 8p . This implies that -1/2-64 describes a particular solution with two arbitrary 
constants. 

For a non-flat brane on the other hand, the dynamical evolution of the braneworlds 
under consideration depends strongly on the form of their bulk matter: although FDSs 
cannot arise when a scalar field fills the bulk, they do become possible once this field is 
replaced with a perfect fluid. In particular, both -1/2-64 and -1/263, with r < — 2 offer the 
possibility of a collapse-brane II singularity, while the balance 7 £>2 with 7 < —1/2 implies 
the existence of a FDS around which a —>■ and a' —>■ a, i.e., a and a' behave exactly as 
in the case of a collapse-brane II singularity but the energy leak into the extra dimension 
is monitored each time by the 7 parameter that defines the type of fluid (or the dominant 
exponent). However, for a non-flat brane there still exists the possibility of escaping the 
FDSs and this is suggested by the balances 7 £>2 with 7 > —1/2 and -1/263 with r > —2. 
(Recall that all possible balances in the case of a non-flat brane correspond to particular 
solutions with two arbitrary constants.) 

It would be interesting to further investigate if the FDSs found here continue to emerge 
when we consider the case in which the scalar field coexists with the perfect fluid in the 
bulk. The analysis of this more involved case may also shed light to the factors that control 
how these two bulk matter components compete on approach to the singularity or even 
predict new types of FDSs that might then become feasible. 



Chapter 6 



Conclusions and future work 



Two main themes throughout this thesis were completeness and singularities. In particular, 
in Chapter 2 we stated and proved the general completeness theorem of [72] for the case 
of an isotropic spacetime. This theorem indicates the underlying reasons that lead to 
completeness in the category of isotropic models, some examples of which we analysed 
later on in that Chapter. 

We continued in Chapter 3 by formulating necessary conditions for the appearance of 
finite-time singularities in isotropic universes based entirely on the behaviour of the Hubble 
parameter. These conditions provided us with a first classification of such singularities, 
a scheme which was then further refined and expanded with the use of the Bel-Robinson 
energy. This energy represents the energy of the gravitational field projected in a sense to 
a slice in spacetime and it has been previously used for the proof of global-in-time results 
(see |206j . |207| . |208j ). Combining the quantities that constitute the Bel Robinson energy 
with the pressure and density of the matter fields through the field equations, we were 
thus led to a clear picture of how the matter fields influence the nature of the singularities. 
The new classification scheme is based on the behaviour of the Hubble parameter, the 
Bel-Robinson energy and the scale factor. As it was shown later on in Chapter 4, an 
examination of models from the recent literature illustrates that this classification not only 
includes all known types of singularities but it also predicts new ones. We believe that 
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this classification is complete in the isotropic category in the sense that it depicts all the 
different qualities in which singularities may differ. Chapter 3 ended with the formulation 
of a relation that holds between the existence of a closed trapped surface and the divergence 
of the Bel-Robinson energy; this relation provided us with the information needed to decide 
on the fate of isotropic cosmologies and we illustrated this important result with various 
examples. 

In Chapter 4, we analysed under the prism of our classification a wide variety of 
isotropic cosmologies that accommodate during their existence one or even more of the 
four highlighted singularity types of our classification scheme, namely a collapse, big- 
rip, sudden, or a sudden-rip singularity. There are many more types predicted by our 
classification and the search for adequate physical circumstances giving rise to these types 
of singularities is a project of future work. 

Incorporating these new singularities in the evolution of particular cosmologies will 
lead to a precise understanding of the role that these singularities may play in determining 
the overall destiny of such cosmologies. Recall from Section 4.3 that it was only after 
observations indicated that dark energy may be described by the phantom equation of 
state that the possibility of a big-rip singularity was considered seriously for the first time. 
This led to the discovery that our universe may face a violent end on approach to which 
every bound structure such as galaxies, stars, planetary systems or even atoms would be 
ripped apart. 

A similar situation may arise in the consideration of models accommodating the not- 
so-far studied singularity types of our classification: in the case that some model comes in 
agreement with observations we may be able to make predictions about the future evolution 
of the observed universe. 

Another project of future work is to extend the present classification by first considering 
anisotropic and then inhomogeneous spacetimes within the framework of general relativity; 
we have already discussed in Chapter 3 (see Section 3.5) the way that the classification 
scheme presented in this work is expected to be enriched with singularities that are by 
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necessity impossible to occur in isotropic models. 

A step towards another direction is to examine how our classification does change when 
we abandon general relativity and consider other metric theories of gravity such as f(R) 
theories or scalar-tensor gravity. With this examination we will be able to recognise the 
way that these different theories of gravity affect the number and types of singularities 
that are possible to occur in a given class of spacetimes. 

The main body of this thesis ended with Chapter 5 at which point we took a turn from 
the classical case of a four-dimensional spacetime to explore the character of singularities 
that arise in the different cosmological landscape of braneworlds. Our research focused on 
the asymptotic analysis of singularities that may arise during the dynamical evolution of a 
three-brane embedded in a five-dimensional bulk space that is inhabited by a scalar field or 
a perfect fluid. Our analysis, which was performed with the method of asymptotic splittings 
of [76], shows that there are three main possible types of finite-distance singularities which 
we call collapse-brane I, II and big rip-brane singularities. A common characteristic of 
these singularities is that they act as a reservoir through which all the energy density of 
the bulk-matter components may be emptied. 

This phenomenon of an infinite leak of energy into the extra dimension was previously 
used in [74J as an explanation to the cosmological constant problem. Their model consisted 
of a flat three-brane embedded in a five dimensional bulk with a scalar field and exhibited 
a collapse-brane I singularity that successfully cut off the extra dimension and swallowed 
the vacuum energy, offering thus an explanation for the observed four- dimensional flatness 
and the vanishing value of the cosmological constant. We showed that this singularity is 
removed to an infinite distance from the brane when instead of a flat brane we considered 
a curved one while by altering the matter in the bulk we allowed for the new distinct types 
of finite-distance singularities mentioned above to emerge. 

Future work in the field of braneworld singularities includes the construction of new 
versions of our model with the consideration of other forms of bulk matter such as the 
combination of a scalar field and a perfect fluid, or a fluid with unconnected pressure and 
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density. Studying different forms of bulk matter may lead to the discovery of new finite- 
distance singularities which will help us to draw conclusions about the possible variety of 
singularities that these braneworlds can accommodate and hence help to determine their 
true range of validity. 



Appendix A 

Notation and conventions 



In this thesis we use the signature (—,+,+,+) for the metric tensor g ab . The Riemann 
tensor is defined by 

R a bcd = d c r a bd - d d ri + r a cf r( d - r^r£, (A.o.i) 

where d a stands for the partial derivative d/dx a and 

= \g ac {d h g cd + d d g cb - d c g bd ). (A.0.2) 

are the Christoffel symbols. The Ricci tensor is defined by contracting the first and the 
third index of the Riemann tensor, i.e., 

R ab = R c acb - (A.0.3) 

The Einstein tensor is given by 

G ab = R ab - ^g ab R, (A.0.4) 

where R is the scalar curvature, R = g ab R ab . 
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In Chapters 2,3,4 we use units so that 8ttG — c — 1. The Einstein equations then read 

G ah = +T ab (A.0.5) 

where T a t, is the energy-momentum tensor. 

The symbol V a stands for the covariant derivative. 

We also use the following abbreviations: RW for Robertson- Walker, FRW for Friedmann- 
Robertson- Walker, CTS for closed trapped surface, FDS for finite-distance singularity and 
IDS for infinite-distance singularity. 



Appendix B 

Some submanifold geometry formulae 



In this appendix we review some useful formulae from submanifold geometry. In Section 
B.l, we calculate the extrinsic curvature of the spacelike slices of an arbitrary globally 
hyperbolic spacetime. In Subsection B.l.l, we show that the extrinsic curvature of an RW 
space is proportional to the Hubble rate. In Section B.2, we derive the basic Gauss-Codazzi 
relations that associate the Riemann and the Ricci tensor of the spacelike slices with the 
Riemann and the Ricci tensor of the spacetime. 

B.l Extrinsic curvature 

Consider a globally hyperbolic spacetime (V, g), that is a spacetime which can be foliated 
by Cauchy surfaces E ( eEx {t} parametrised by a global time function t, and a smooth 
congruence of timelike geodesies in (V, g) orthogonal to the spacelike hypersurfaces S t . If 
n a is the unit tangent vector field (n a n a = —1) to these geodesies, then n a is also the unit 
normal vector field to the hypersurfaces E t . The spacetime metric g ab induces the spatial 
metric h ab on £ t given by [77] 

hab = g a b + n a n b . (B.l.l) 

It then follows that h ab satisfies 
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h ab n a = 0, h ab h b c = h ac . (B.1.2) 

To understand the role that h a b plays consider a spacetime vector w a . We can split w a 
into components tangent and normal to E t so that 

w a = w ± n a + wl (B.1.3) 

where w^n a = 0. If w a is a tangential vector, then w a = and w± = 0, or equivalently 

w a = h a b w b . (B.1.4) 

This condition can be generalised to apply to an arbitrary tensor in the following way: we 
say that the tensor T£^"£ k is tangential to S 4 if 

= Kl ■ ■ ■ KtK ■ ■ ■ < ; T-;.t • (b.i.5) 

This means that h% is the projection operator onto the tangential plane perpendicular to 
n a . 

The extrinsic curvature K a b of St is defined as 

K ab = V a n b . (B.1.6) 

Therefore the extrinsic curvature describes the bending of E 4 in the spacetime. To see this 
more clearly, consider the orthogonal deviation vector £ a representing the spatial infinites- 
imal displacement between two neighboring geodesies of the congruence. Then £ n £ a = 0, 
which means that = K^ b . The last equality implies that the extrinsic 

curvature measures the failure of n a to be parallelly transported along a spacelike geodesic 
with tangent vector £ a . In the same time, the first equality says that this failure is equal 
to the failure of £ a to be parallelly transported along the timelike geodesic with tangent 
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n a . This means that an observer on one of the geodesies would see a nearby geodesic being 
stretched and twisted [209J. 

Note that K a b is purely spatial because K a bn b = 0. Since we have considered a 
hypersurface-orthogonal congruence of geodesies, it can be proved (cf. [77], p. 217) that 
the rotation u a b, defined by uj a b = \iSbn a — V a rib) = K^a], is zero and hence the extrinsic 
curvature tensor is symmetric K a b = Kb a - We can therefore write 

2K ab = V a n b + V b n a . (B.l. 7) 

The right hand side of the above equation equals the Lie derivative of g a b with respect to 
n a , hence 

11 11 

K ab = -£n9ab = -£ n {h a b - n a n b ) = -C n h ab - -C n (n a n b ). (B.l. 8) 

But 

C n (n a )n h = (n c V c n a + n c V a n c )n b = (^g a dn c V c n d + ^V a (n c n c )^ n b = 0, (B.l. 9) 

because of the geodesic equation n c V c n d = 0, and the fact that n c n c — — 1. Using this and 
the Leibnitz rule in Eq. (IB. 1.81) . we find the following formula for the extrinsic curvature: 

K ab = \C-nhab- (B.1.10) 

Since a Gaussian coordinate system is adapted to n a , we finally get 

K ab =-^. (B.l.ll) 

ab 2 dt y 1 



This relation states that the extrinsic curvature is the rate of change of the spatial metric 
h a b along the congruence. 
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B.l.l Extrinsic curvature for a RW space 



Consider now a RW metric, i.e., one with 



ds 2 = -dt 2 + a 2 (t) (dr 2 + f 2 (r)(d6 2 + sin 2 Odcj) 2 )) 



(B.1.12) 



where 



fir) = { 



r, k = 

sinr, k — 1 
sinhr, k = — 1. 



Then by the use of Eq. ( IB. 1.111) . we find 



(B.1.13) 



K\i = a(t)a(t), 

K 22 = f 2 (r)a(t)d(t), 

K 33 = f 2 (r) sin 2 9 a(t)a(t). 

Now, since the spatial norm of a 2-covariant tensor is given by 



\X\ 2 = g i *g kl X ik X jl i,j,k,l = 1,2,3, 



(B.1.14) 



we find that the norm of the extrinsic curvature is given by 

\K\ 2 = g^f s K ir K js = (g^ 2 K 2 n + (g 22 ) 2 K 2 22 + (g^) 2 K 2 3 = = 3H 2 (t), (B.1.15) 

where H(t) is the Hubble expansion rate, H = a/a. 



B.2 Gauss-Codazzi formalism 



As we saw in the previous Section, the metric g a b induces a Riemannian metric h a b on each 
hypersurface S 4 . We know that there exists a unique derivative operator associated with 
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this metric, which we denote by D c . It can be proved that D c is given by [77] 

DcT^X" = Ki . . . K]h{v f Ttl;;i k ■ (B.2.1) 

In turn D c determines a curvature tensor ^R d bc on S t . This tensor is defined in analogy 
with the four dimensional case by the relation 

D a D b u c - D b D a u c =^ R d abcUd . (B.2.2) 

Our goal is to find the relation between the Riemann tensor of the hypersurface ^R^ 
and the Riemann tensor of the whole spacetime R d bc (Eq. IB.2.121 below). We start by 
calculating D a D h u c : 



D a D b u c = D a (hXV d u e )=hlh k XV f (h e k h d g V d u e ) 

= hih k c hth d y a (v f v d u e ) + h d g hih k c h{iy d u e )v f K 

+ hlh k c hiht(V d u e )(V f h d g ). (B.2.3) 



But by the use of 



h{h k c V f K = h{h k V f {gt + n k n e ) 

= h{h k c {V f n k )rf + h{h k c n k V s rf + h{h k V f g e k 

= h k K ak n e = K ac n e , (B.2.4) 



Eq. (IB.2.3P becomes, 



D a D b u c = h 9 b h d h k hlh{{V s V d u e ) + K ac n e h d (V d u e ) + h e c K ab n d (V d u e ). (B.2.5) 



Now since 



h d iy d u e )n e = hy d (u e n e ) - h d (V d n e )u e = -K e b u e , 



(B.2.6) 
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we find for D a D b u c , 

D a D b u c = hlhfhy f V d u e + h e c K ab n d V d u e - K e b u e K ac . (B.2.7) 
Similarly, we get an analogous expression for D b D a u c , 

D b D a u c = h d b h{hy d V f u e + h e c K ba n d V d u e - K e a u e K bc . (B.2.8) 

Subtracting Eqs. (IB.2.3D and flB.2.81) we have that 

D a D b u c - D b D a u c = -K ac K^u e + K bc K e a u e + h f a h d b h e c {VfV d u e -V d V f u e ) 
= -K ac K e b h d u d + K bc K e a h d u d + h{h d h e c R) de u 3 
= -K ac K d u d + K bc K d u d + h{h d hlh d R j fde u d . (B.2.9) 

Then from Eq. (IB. 2. 91) and Eq. (1B.2.2j) follows the first Gauss-Codazzi relation: 



(3) ^L = h{h d Khp fde - K ac K d + K bc K d . (B.2.10) 

In Gaussian normal coordinates it holds that 

H = hf x h Xa = 5l (B.2.11) 

so that Eq. (|B.2.10[) becomes 

{3) R d abc = R d abc ~ K ac K d + K ab K d , (B.2.12) 



where we have used the fact that K bc K d = gi c g la K ab K d = K ab K t 



Now we can use the Eq. (IB. 2. 121) in order to deduce the analogous relation that holds 
between the Ricci tensor of the hypersurface ^R a b and the Ricci tensor of the spacetime, 
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R ab . For that purpose, we start by calculating the component R™ nb (cf. |210j . pp. 163-164): 



b t.a 

m 



R d arb n d n r h b q h a m = (V a V r n b - V r V a n b )n r h b q h a m = (V a K rb )n r h b q h a m - (V r K ab )n r h b q h 

= V a {K rb n r h b q h a m ) - K rb (V a n r )h b q h a m - K rb n r V a (h b q h a J 

- (£ n K ab - K rb V a n r - K ar V b n r )h b q h a m 

= (K ar K r b - £ n K ab )h b q h a m , (B.2.13) 



so that 



R d arb n d n r = K ar K r b - £ n K ab . (B.2.14) 



Introducing Gaussian-normal coordinates and contracting the first and the third index, 
we have that 

R n anb = ~K ar Kl + ^. (B.2.15) 



The Ricci tensor is then given by 



dK ab 



Rab = Kmb + Knb = (3) Rat + K ab K - 2K an K£ + (B.2.16) 

on 

We can work in a similar way for the case of a timelike hypersurface with normal vector 
n a (now n a n a =l). Both cases of a spacelike and timelike normal vector can be treated at 
the same time throughout the previous calculations by using the symbol e(n) = n a n a = ±1. 
In particular, with the use of e, the key equations flBTTD . flB.2.12|t . (IB.2.14D and flB.2.16|> 
can be written in the following form: 

Kb = 9ab - e(n)n a n b , (B.2.17) 

^R d abc = R d abc + e{n){K ac K d b -K ab K d c ), (B.2.18) 

R n anb e{n) = K ar K r b - £ n K ab , (B.2.19) 

Rab = {3) R ab + e(n)^-K ab K + 2K an KZ-^y (B.2.20) 
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Appendix C 
Bel-Robinson energy 



This appendix deals with the notion of the Bel- Robinson energy. In Section C.l, we give 
the definition of the Bel-Robinson energy for the general form of a sliced space and in 
Subsection C.l.l, we calculate the Bel- Robinson energy of a RW space. In the end of 
the Subsection C.l.l, we show that the null energy condition in closed or flat accelerating 
universes translates into an inequality satisfied by the electric parts of the Bel-Robinson 
energy. 

C.l Definition 

The Bel-Robinson energy and related quantities play an important role in proving global 
properties of spacetimes. For example, they were used in [208J, [207J to prove global 
existence results in the case of asymptotically flat and cosmological spacetimes respectively. 
In this thesis, the Bel-Robinson energy is used for the purposes of classification of spacetime 
singularities. 

The Bel-Robinson energy is a kind of energy of the gravitational field projected in a 
sense to a slice in spacetime. Consider the sliced spacetime with metric 

( 3+1 ^ = -N 2 {6 ) 2 + g i3 Q l Q\ 6° = dt, l = dx l + fidt, (C.l.l) 
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where N = N(t, x l ) is the lapse function and /3 l (t, x- 7 ) is the shift function. The 2-covariant 
spatial electric and magnetic tensors are defined as follows 

E ij = R iOj> (C.1.2) 

Dij = -VihkVjimR hklm , (C.1.3) 
Hn = l -N- l r] ihk Rl^ (G.1.4) 
Ba = ^N-\ hk R^, (C.1.5) 

where rjijk is the volume element of the space metric g% = gijdx l dxK These four time- 
dependent space tensors comprise what is called a Bianchi field, (E, H, D, B), a very im- 
portant frame field used to prove global in time results, cf. [01]. The Bel-Robinson energy 
at time t is given by 

B(t) = ~ / (\E\ 2 + \D\ 2 + \B\ 2 + \H\ 2 ) dfi 9t , (C.1.6) 
where by \X\ 2 = g^ g kl Xi k Xji we denote the spatial norm of the 2-covariant tensor X. 

C.l.l Bel- Robinson energy of an isotropic universe 

Consider the RW metric (I2.1.2p . In order to compute the Bel- Robinson energy at time t, 
Bit), we first calculate the two covariant spatial tensors Ey and using relations (1C.1.2j) 
and (lC.1.3p . These relations give 

En — oa, 

E 2 2 = fad, 

E 33 = f 2 sin 2 9 ad, 



while 
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{ f" 

2 / a2 J 



D 33 = /Wflfd 2 -^ 



where a prime denotes the differentiation of / = f(r) with respect to r, while a dot denotes 
the differentiation of a = a(£) with respect to i. The magnetic tensors f/y, Bji make no 
contribution to the Bel-Robinson energy since for a RW metric we have that = for 
h,k,j = 1,2,3. 

Next, we make use of Eq. (IB. 1. 141) in order to calculate the spatial norms of the two 
electric tensors. We find 

\E? = (g n ) 2 E 2 u + {g 22 ?El 2 + {g^fEt, = 3^, (C.1.7) 

and 

w = (j 11 ) 2 ^ + (/r^ + (g 33 r^s = (^f + + 2 g - . (c.1.8) 

It then follows from the form of f(r) (cf. Eq. (12.1.31) ) that 

1 - f 2 f" 

where /c is the constant curvature of the spatial slice. We can therefore write \D\ in the 
form 

|i>l» = 3 (*+*)'. (C.1.10) 

The Bel-Robinson energy at time t is then given by 
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e(() 4(|^ + | D | 2) ^g + (^ + *V), (0111) 



where C is the constant volume of (or in, in the case of a non-compact space) the three 
dimensional slice at time t. 

Finally, we have the following proposition that translates the null energy condition 
into an inequality satisfied by the electric parts of the Bel-Robinson energy. 

Proposition C.l.l In an accelerating flat or closed universe the null energy condition is 
equivalent to the inequality \E\ < \D\. 

PROOF. The null energy condition demands that 

R ab k a k b > 0, for every null vector k a . (C.1.12) 

In an isotropic universe with metric described by Eqs. (12.1.2ft and (12.1. 3D . the above 
condition translates to 

R ab k a k b = R^k 1 ) 2 + R 22 (k 2 ) 2 + R 33 (k 3 ) 2 + R u (k 4 ) 2 

nil \ 

;2 nJ , „■■ \ /;„1\2 , f-t jt/2 r nil , f2/ -2 , "\\/;„2\2 



2d 2 - 2y + adj (A; ) + (1 - f' 2 - ff" + f 2 (2d 2 + ad))(k 2 ) 2 + 

+ sin 2 6(1 - f 2 - ff" + f 2 {2d 2 + ad))(k 3 ) 2 - 3-(k 4 ) 2 . (C.1.13) 

a 

Keeping in mind that 1 — f' 2 = —ff", we can write Eq. (1C. 1.131) as follows 

R ab k a k» = (2^-+'-- 2-0 [a 2 ({k 1 ) 2 + f 2 ((k 2 ) 2 + sin 2 9(k 3 ) 2 )) - (fc 4 ) 2 ] + 



( -3- + 1% + - ~ 2-£) (k 4 ) 2 = (-2* +2%-2^-\ [k iVI 



a? a fa? 
a d 2 a 

!- + 2 h - - ,.. , - 

a a 2 a fa 2 J \ a a 2 fa 

Therefore if the null energy condition, Eq. (10. 1 . 12|) . holds then 



a a k , . 

- + ^ + ^>0, (C.1.14) 



a a 2 a 2 
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where we have used the fact that — /"// = k, with k being the constant curvature of the 
spatial slices. Since we have assumed a k > accelerating universe, the electric parts are 
given by 

LB = v^- 
a 

hence Eq. flC. 1.14ft implies that \E\ < \D\. m 
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Appendix D 



Energy-momentum tensor of a scalar 
field 



In order to define the energy-momentum tensor for a scalar field we start from the action 



;d.o.i) 



where the Lagrangian density is 



C-4, = -^(V c 0V c + 2V(<f>)) = ~{g dc V d <pV c <P + 2V (((>)). 



(D.0.2) 



The energy-momentum tensor for a matter field is defined to be, cf. |77J, p. 455, 



T, 



2 S(£ m atter \J ~ 9) 

f-9 &9uv 



;d.o.3) 



for the Lagrangian (ID. 0.21) . we find that 



T, 



-9 



±= ((V c 0V c + 2V(<f>))(~V=99^) + V^VAyfH 



(D.0.4) 
(D.0.5) 



115 



116 APPENDIX D. ENERGY-MOMENTUM TENSOR OF A SCALAR FIELD 

Thus the energy-momentum tensor of a scalar field is 

Tfa, = V M 0V,0 - ^(V c 0V c + 2V(<j>)), (D.0.6) 



or equivalently, 



= V M 0V^ - \g^{g^ d (jN c (j) + 2V (</>)). (D.0.7) 



Comparing this with the energy-momentum tensor of a perfect fluid, 

T^v = {p + p)u^u u + pg^, (D.0.8) 
where is the unit timelike vector (u^u^ = —1), it follows that 

P = C 4> = -i(^ c V d 0V c + 2V ((/>)), (D.0.9) 

and 

(jjl + p)u^u v = V^V„0. (D.0.10) 
Contracting the above equation, we find (99] 

-{n + p) = V>V„0 = -2L^ - 2V(4>), (D.0.11) 

and substituting for the pressure from Eq. flD.0.9j) leads to 



/i = C <P + 2V((P). (D.0.12) 



In a homogeneous universe, a scalar field is a function of time only and hence Eqs. flD.0.9j) 
flD.0.12j) give the standard result, 



P = \-V{4>), (D.0.13) 
/i = ^+^(0). (D.0.14) 
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The equation of motion of a scalar field can be found either by using the conservation law 
of energy, i.e., V u T liU = 0, or by variation of the action with respect to 0. We use the 
second approach. We have that: 

8Cs dCs ( d£ 



(^v)) = 2y ' (0) " 2VcVC0 ' (D '°' 15) 



5<j> d<j> VWA, 

and so 5 C^/ 5(f) = gives directly 

□0 - V\(f>) = 0, (D.0.16) 
where U<f> = (-g)-^V d ((- g y/ 2 g dc V c( f>). 
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Appendix E 

The method of asymptotic splittings 



Consider a system of n first order ordinary differential equations 



x = f(x), 



(E.0.1) 



where x = (xi, . . . ,x n ) G M n and f(x) = (fi(x), . . . , / n (x)). The general solution of the 
above system contains n arbitrary constants and hence describes the behaviour of the 
system for arbitrary initial data. A particular solution on the other hand contains less 
than n arbitrary constants and describes the behaviour of the system for a subset of initial 
data. 

We say that a general or particular solution of the above dynamical system exhibits a 
finite-time singularity if there exists ai s 6l and a xq G W 1 such that 



where || . || is any L p norm. 

The purpose of singularity analysis (cf. [75] . |211j . pp. 129-155) is to build series 
expansions of solutions around the presumed position of a singularity at t s in order to 
study the different asymptotic behaviours of the solutions of the system (IE. 0.11) . 

In particular, we look for series expansions of solutions that take the form of a Puiseux 




(E.0.2) 
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series, 



x = r p (a + Y,™ lCi T i/s ) , (E.0.3) 

where r = t — t s , p E Q n , sGN. 

The method of asymptotic splittings for the system (lE.O.ip is realised in the following 
stepsj: 

• First, we find all the possible weight-homogeneous decompositions of the vector field 
f by splitting it into components f^: 

f = f(°)+f( 1 ) + ... + f(*) j (E.0.4) 
with each of these components being weight homogeneous, that is to say 

f(J)( ar P) = r p+i(9«-i) f 0-)( a ) j = o, . . . , k, (E.0.5) 



where aGR" and are the positive non-dominant exponents that are defined by (1E.0.7j) 
below. 

• We then substitute the forms x = ar p in the system x = f(°)(x) in order to find 
all possible dominant balances, i.e., sets of the form {a, p}. The order of each balance is 
defined by the number of the non-zero components of a. 

• For each of these balances we check the validity of the following dominance condition 

lim = 0, (E.0.6) 

and define the non-dominant exponents j — 1, . . . , k by the requirement 

E* =1 fCi)( r p 



-p-i 



~r 9<J) . (E.0.7) 



The balances that cannot satisfy the condition (IE.0.6[) are then discarded 



1 Here we will refer briefly to the basic steps of the method of asymptotic splittings. A detailed analysis 
can be found in [76], [212] . 
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• Next, we compute the Kovalevskaya matrix /C defined by 

/C = Df (0) (a) -diagp, (E.0.8) 

where Df(°> (a) is the Jacobian matrix of f ^ evaluated on a. 

• We calculate the spectrum of the /C-matrix, spec(JC), that is the set of its n eigen- 
values, called also the JC-exponents. The arbitrary constants of any particular or general 
solution first appear in those terms in the series (1E.0.3j) whose coefficients have indices 



k = ps, where p is a non-negative /C-exponent and s is the least common multiple of the 
denominator of the set consisting of the non-dominant exponents qv> and of the positive 
/C-exponents (cf. [76], |212j ). The number of non-negative /C-exponents equals therefore 
the number of arbitrary constants that appear in the series expansions of flE.0.31) . There 
is always the —1 exponent that corresponds to an arbitrary constant that is the position 
of the singularity at t s . A dominant balance corresponds thus to a general solution if it 
possesses n — 1 non-negative /C-exponents (the nth arbitrary constant is the position of the 
singularity, t s ). 

• Next, we substitute the Puiseux series: 



Xl = Ef =oCji r Pi+ ^ s , i = l,...,n, (E.0.9) 



in the system flE.O.lj) 



• We find the coefficients Cj by solving the recursion relations 

Kcj - -Cj = Pj-(ci, . . . , Cj-_i) (E.0.10) 

where Pj are polynomials that are read off from the original system. 

• We verify that for every j = ps, with p a positive /C-exponent, the following com- 
patibility conditions hold: 

v T -Pj = 0, (E.0.11) 
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where v is the eigenvector associated with the positive /C-exponent p. 
• We repeat the procedure for each possible decomposition. 

We note here that if the compatibility condition is violated in some eigenvalue then 
the Puiseux series is not valid and instead we have to substitute a ip -series, which is a 
series of the form 

x = r p (a + S^S^c^-r^Klogrp) , (E.0.12) 

where p is the first /C-exponent for which the compatibility condition is violated. The rest 
of the procedure in this case is the same as before. 
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